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Abstract 

A large number of multifaceted quantum transport processes in 
molecular systems and physical nanosystems, such as e.g. nonadia- 
batic electron transfer in proteins, can be treated in terms of quantum 
relaxation processes which couple to one or several fluctuating envi- 
ronments. A thermal equilibrium environment can conveniently be 
modelled by a thermal bath of harmonic oscillators. An archetype sit- 
uation provides a two-state dissipative quantum dynamics, commonly 
known under the label of a spin-boson dynamics. An interesting and 
nontrivial physical situation emerges, however, when the quantum dy- 
namics evolves far away from thermal equilibrium. This occurs, for 
example, when a charge transferring medium possesses nonequilib- 
rium degrees of freedom, or when a strong time-dependent control 
field is applied externally. Accordingly, certain parameters of underly- 
ing quantum subsystem acquire stochastic character. This may occur 
for example for the tunnelling coupling between the donor and ac- 
ceptor states of transferring electron, or for the corresponding energy 
difference between electronic states which assume via the coupling 
to the fluctuating environment an explicit stochastic or deterministic 
time-dependence. Herein, we review the general theoretical frame- 
work which is based on the method of projector operators, yielding 
the quantum master equations for systems that are exposed to strong 
external fields. This allows one to investigate on a common basis the 
influence of nonequilibrium fluctuations and periodic electrical fields 
on those already mentioned dynamics and related quantum transport 
processes. Most importantly, such strong fluctuating fields induce a 
whole variety of nonlinear and nonequilibrium phenomena. A charac- 
teristic feature of such dynamics is the absence of thermal (quantum) 
detailed balance. 
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1 Introduction 

The description and analysis of the dynamics of open quantum systems, 
i.e. quantum systems interacting with a dissipative environment, presents 
a key challenge for nonequilibrium statistical physics. Moreover, this theme 
is also of prominent importance for many applications in physics, physical 
chemistry and physical biology. This can be exemplified by the relaxation 
dynamics occurring in a two-level quantum system that is coupled to the 
vibrational degrees of freedom of an environment. This latter theme gained 
great popularity and is known under the label of " spin-boson dynamics" [1-4] . 
Several apparently different physical problems can formally be unified within 
such a common mathematical description. For example, the relaxation dy- 
namics of a nuclear spin 1/2 in solids, the tunnelling of defects in metals, 
the relaxation dynamics of atoms in optical cavities can all be modelled by 
a (pseudo)-spin-boson dynamics. Another important and relevant situation 
refers to donor-acceptor electron transfer reactions in various molecular struc- 
tures [5-9]. For spatially extended quasi-periodic molecular structures like 
those formed by protein a-helices [10-13], or DNA's [14-18] many quantum 
states are generally required to describe charge transfer processes. Here, 
a multi-state tunnelling problem naturally emerges with the tight-binding 
model serving as one of simplest theoretical frameworks. 

The primary problem is to describe the influence of the environmental 
degrees of freedom on the quantum dynamics of interest. Many different 
approaches have been developed to tackle this challenge. The fundamental 
methodology consists in separating the total system into the two (or more) 
mutually interacting parts: the dynamical subsystem with a small number 
of relevant degrees of freedom and a thermal bath represented by a huge 
number of microscopic degrees of freedom which are at thermal equilibrium. 
A most general quantum-mechanical description is provided by the density 
operator of the whole, combined system which depends both on the variables 
of the considered dynamical subsystem (relevant variables) and the variables 
of the thermal bath (irrelevant variables). The dynamical behaviour of a 
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small quantum subsystem presents the focus of interest with the thermally 
equilibrated bath degrees of freedom serving as a source of randomness for 
the relevant dynamics. This randomness can effectively be eliminated via a 
course-grained description of the system of interest. A corresponding aver- 
aging procedure results in a contracted, reduced dynamics which generally 
entails memory effects, decoherence and dissipation. 

Different approaches have been developed over the years within this gen- 
eral line of reasoning. Within a variety of different approaches, the method 
of path integrals in real time [1,2,19-21] and the projection operator method 
[22-31] provide some of the most frequently used methods. The path-integral 
approach can, however, be technically cumbersome in practical applications 
of interest. The projection operator method appeals to its generality and 
technical elegance. It allows one to obtain formally exact generalised master 
equations (GMEs) for the reduced density matrix in a straightforward way. 
By and large, however, such exact GMEs cannot be analytically elaborated 
further without invoking some sort of a perturbation technique with corre- 
sponding approximations. For example, already the seemingly simple spin- 
boson dynamics cannot be solved analytically exactly. The weak-coupling 
approximation of the system-bath coupling is one of the most useful and 
commonly employed scheme. Moreover, a strong-coupling problem can of- 
ten be mapped onto a (different) weak-coupling problem within a canonically 
transformed basis of the total system. The projection operator method, com- 
bined with appropriate canonical transformations, further improved by use 
of variational approaches, presents a powerful and general method of wide 
acceptance. This well-established methodology is, however, also rather de- 
manding. 

Yet another popular methodology consists in modelling the thermal bath 
influence through a classical stochastic field which acts upon the considered 
dynamical system. Formally, this methodology corresponds to introduction 
of randomly fluctuating time-dependent forces in the Hamiltonian of con- 
sidered quantum system [32-35] and finding subsequently the stochastically 
averaged evolution of the considered system which is governed by a stochas- 
tic Liouville-von-Neumann equation. This methodology is known under the 
label of stochastic Liouville equation (SLE) approach [34-41]. Due to a rea- 
soning that involves the central limit theorem classical random forces with 
Gaussian statistics are most frequently used in this kind of approximate mod- 
elling. The Gaussian white noise serves here as a simplest implementation for 
the corresponding classical stochastic bath. It corresponds to a bath with an 
infinite spectrum of excitations. Such models can be solved exactly in a num- 
ber of specific cases [35,38,40-44]. All the thermal baths possess, however, 
finite energy spectra. This circumstance gives rise to temporal autocorrela- 
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tions in the bath-induced classical stochastic fields. Gaussian Markov noise 
with the exponentially decaying temporal autocorrelations presents one of 
the simplest models of such more realistic, coloured noise [45]. Yet, even in 
the simplest case of a two-state tunnelling system this model cannot be solved 
exactly except for some limiting cases (see, e.g., in [46] for the Landau-Zener 
model with a stochastic modulation). Generically one must invoke some ap- 
proximations; e.g., in the case of a weakly coloured Gaussian noise some kind 
of cumulant expansion technique [38,39,47] can be used. 

There exists a different possibility. Continuous state noises can be ap- 
proximated by noise sources with a large number of discrete states (e.g., by a 
discretisation procedure of a continuous diffusion process in a potential). Cer- 
tain Markovian discrete state noises provide then a rather general framework 
for a formally exact stochastic averaging [34,48,49]. Moreover, the two-state 
Markovian noise (also known as dichotomous noise) presents such a simple 
discrete noise source which allows for an exact study of noise driven two- 
level quantum systems [45,50-56]. In addition, the multistate case of exciton 
transfer in molecular aggregates with many quasi-independent noise sources 
modelled by independent two-state Markovian noises presents another ana- 
lytically tractable case, in the sense that it can be reduced to the solution of a 
system of linear differential equations with constant coefficients for averaged 
dynamics [57]. The discussed dichotomous noise can model a quasi-spin 1/2 
stochastic bath variable. In the case of electron transfer in molecular systems 
such a quasi-spin stochastic variable can simulate, for example, the bistable 
fluctuations of a charged molecular group nearby the donor, or acceptor site, 
or the conformational fluctuations of a bistable molecular bridge. 

A well-known drawback of the SLE approach consists, however, in the 
asymptotic equal-population of the energy levels of quantum system which 
occurs for arbitrary energy differences [38,40,41]. This means that the SLE 
approach corresponds formally to an infinite bath temperature. At least, 
the thermal energy ksT should thus be larger than the characteristic en- 
ergy scale of the quantum system, e.g., larger than the energy width of the 
corresponding excitonic band. This corresponds to a high-temperature ap- 
proximation [40,41,52,57]. The reason for this intrinsic restriction is that 
the stochastic field unidirectionally drives the quantum system without being 
modified by the system's feedback (no back reaction). This drawback within 
the SLE approach requires some ad hoc corrections to enforce the correct 
thermal equilibrium [42,58,61]. Nevertheless, the SLE approach can yield 
a very useful tool, notably in the nuclear magnetic resonance (NMR) the- 
ory [40,59], the theory of exciton transfer in molecular aggregates [40] and 
within the theory of single-molecular spectroscopy [60]. 

Combined scenarios have been used in several works [61-63]. Initially, 
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those were aimed to model the influence of relaxation processes in the ther- 
mal bath [61], or to account for non-Gaussian large- amplitude fluctuations 
of molecular charged groups [62]. However, it was soon recognised that the 
addition of classical noise to a dissipative quantum dynamics generally vio- 
lates the detailed balance symmetry at the environmental temperature [64]. 
Therefore, the stochastic field in these approaches correspond physically to 
a nonequilibrium noise influence. It has been shown theoretically that such 
a nonequilibrium non-Gaussian (e.g., two-state) noise can regulate the quan- 
tum transition rates by several orders of magnitude [62,64-66]. Moreover, it 
may pump energy into the quantum system. This in turn gives rise to various 
interesting nonlinear nonequilibrium phenomena such as a noise-induced en- 
hancement of thermally assisted quantum tunnelling [67], an inversion of pop- 
ulation in discrete quantum dissipative systems [68] , a noise- induced absolute 
negative mobility for quantum transport [69], or also a fluctuation- induced 
transport of quantum particles within a tight-binding description [70], to 
name the most prominent ones. From a thermodynamical perspective these 
nonequilibrium effects are due to a virtual presence of two heat baths of 
different nature: a first one assuming the temperature of the environment 
T (modelled by a thermal bath of harmonic oscillators that are bi-linearly 
coupled to the relevant system), another one possessing a virtually infinite 
temperature T a = oo (stochastic bath). In this intuitive picture, a nonequi- 
librium stochastic field is expected to heat the quantum-mechanical degrees 
of freedom, causing various, surprising nonequilibrium phenomena. 

The study of the dynamics of such quantum dissipative systems which 
are driven far from thermal equilibrium by nonequilibrium fluctuations is 
the focus of present work. The situation here is similar in spirit to one in 
the recently emerged field of (classical) Brownian motors [71-76], see, e.g., 
in [77-81] for surveys and further references. 

In this review, we present a general outline with many important exam- 
ples given of the following methodology: The nonequilibrium stochastic field 
is represented by an external time-varying classical field in the Hamiltonian 
of the quantum system. This field is treated without invoking any further 
approximation, until it becomes necessary to do so. In doing so, a formally 
exact generalised master equation is obtained which includes the external 
field both in the dynamical part and in the dissipative kernel of the GME 
exactly. Subsequently, the dissipative kernel is expanded to the lowest order, 
i.e. the second order in the system-bath coupling. [In a properly canoni- 
cally transformed basis this scheme allows one to study the opposite limit 
of strong dissipation/weak tunnelling as well.] The overall procedure results 
in an approximate, generalised master equation for the reduced density ma- 
trix of the considered, relevant quantum system. We recall that within this 
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methodology the external field is not only included exactly in the dynamical 
part, but it modifies as well the dissipative kernel in a very profound man- 
ner. In particular, the dissipative kernel becomes a retarded functional of 
the driving field. Thereby, the field influence on quantum dynamics is taken 
rigorously into account within the given order of the system-bath coupling. 
Such a corresponding modification of the dissipative kernel becomes crucial 
for strong driving fields. 

This so tailored approach allows one to describe stochastic and time- 
periodic fields on equal grounds. The influence of a time-periodic driving 
on the dissipative quantum dynamics has been investigated in Refs. [68, 
82,83,85-87]. Other related work has been done at the same time and 
in parallel by several other research groups elaborating similar [88-90], or 
different [84,91-97] approaches which reconcile within some approximations 
[3,4,87,98,99]. 

The difference between the stochastic fields and the periodic fields enters 
within our approach on the level of averaging of the corresponding field-driven 
generalised master equations. In this context, one must refer to some further 
approximations, which generally are based on the separation of time scales 
involving the external driving, the contracted quantum dynamics, and the 
decay of dissipative kernels in the generalised master equations. Remarkably, 
in the case of dichotomous fluctuations this averaging can be done exactly 
without further approximations [67,100,101]. The spin-boson model driven 
by such dichotomous Markovian fluctuations presents one instance of general 
interest which has been studied in detail in Refs. [67,68,85,86,101]. Two 
other important situations, where an exact averaging is feasible for a broad 
class of stochastic and periodic processes, are given by an infinitely extended 
tight-binding model. These are: (i) the case of coherent tunnelling in the 
absence of dissipation and (ii) the regime of incoherent tunnelling (strong 
dissipation) when the tunnelling is weak (high tunnelling barriers). Some 
explicit pertinent examples are discussed in Refs. [70,102-104]. 

This review is organised as follows, see in Fig. ^ The study of a quantum 
dynamics subjected to non-Markovian stochastic fields that are modelled by 
discrete state processes of the renewal type is presented in Section 2. Therein, 
a formally exact averaging of the quantum evolution over the stationary re- 
alizations of stochastic fields is given. The general results are illustrated 
by a new Laplace-transformed exact solution of averaged two-level quantum 
dynamics driven by a symmetric non-Markovian two-state field. The prior 
results for a quantum two-level dynamics driven by a dichotomous Markovian 
field are reproduced as a particular limiting case. This section contains also 
the results for the fluctuating Kubo oscillator. These are used in the follow- 
ing. This section also contains a short primer into the projection operator 
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Figure 1: A flow diagram depicting the various topics and their mutual 
interrelations which are covered by this review. 
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formalism. Two-state quantum dynamics in strong periodic fields is consid- 
ered in Sec. 3. The Section 4 outlines the general formalism of dissipative 
quantum dynamics in strong, time-varying fields within the reduced den- 
sity matrix approach. The corresponding weak-coupling generalised master 
equations and the generalised Redfield equations are presented there. These 
equations serve as a basis for subsequent applications and analysis of stochas- 
tic field-induced phenomena. Section 5 contains a simple implementation of 
our general approach which manifests the origin and basic features of strongly 
nonequilibrium phenomena described in the subsequent Sections for more re- 
alistic models. The stochastically and periodically driven spin-boson model 
is discussed in Section 6, including quantum stochastic resonance features. 
Section 7 is devoted to the phenomenon of noise-induced, absolute negative 
mobility in quantum transport and to the analysis of dissipative quantum 
rectifiers. Concluding remarks are drown in the Section 8. 



2 Quantum dynamics in stochastic fields 
2.1 Stochastic Liouville equation 

To begin, let us consider an arbitrary quantum system with a Hamilton 
operator if [£(£)] which depends on a classical, noisy parameter £(£). This 
stochastic process £(t) can take on either continuous or discrete values. Ac- 
cordingly, the Hamiltonian H acquires randomly in time different operator 
values /f [£(£)] which generally do not commute, i.e., [if if [£(i')]] ^ 0. 

The posed problem is to average the corresponding quantum dynamics 
in the Liouville space, which is characterised by the Liouville-von-Neumann 
equation 

j t p{t) = -icmw) , (i) 

for the density operator p(t) over the realizations of noise £(£). £[£(£)] 
in Eq. ((TJ) stands for the quantum Liouville superoperator, £[£(£)](•) = 
}r[.H"[£(t)], (■)]. In other words, the objective is to evaluate the noise-averaged 
propagator 

fto+t 

(S(t + t,t )) = (TexpH / mr)}dr]), (2) 
where T denotes the time-ordering operator. 
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t tj t 2 



Figure 2: Typical trajectory of the considered process [105]. 



2.2 Non-Markovian vs. Markovian discrete state fluc- 
tuations 

We specify this task for a discrete state noise with iV states & (cf. Fig. |2J). 
The noise is generally assumed to be a non-Markovian renewal process which 
is fully characterised by the set of transition probability densities for 
making random transitions within the time interval [r, r + dr\ from the state 
j to the state %. These probability densities are obviously positive and do 
obey the normalisation conditions 

M*r)dT=l, (3) 

i=i Jo 

for all j = 1,2,..., JV. 

The subsequent residence time-intervals between jumps are assumed to 
be mutually uncorrelated. The residence time distribution (RTD) i[ ) j{ T ) i n 
the state j reads 

i 

The survival probability $j(t) of the state j follows then as 

^{r)dr. (5) 

This constitutes the general description for continuous time random walk 
(CTRW) theory [106-109]. 

Several descriptions used for such simplest non-Markovian processes of 
the renewal type are worth mentioning. The approach in Ref. [110] with the 
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time-dependent aging rates ky(t) for the transitions from state j to state i 
corresponds to a particular choice 



ipij(r) := %(r)exp[-2_^ / kij(t)dt\. (6) 



The Markovian case corresponds to time-independent transition rates kij(r) = 
const. Any deviation of ipij{ T ) from the corresponding strictly exponential 
form which can be accounted for by introducing a time-dependence of the 
transition rates kij(r) amounts to a non-Markovian behaviour 1 . Furthermore, 
the survival probability $j(t) in the state j is determined by 

N 

$,(r)=exp[-^ / k l3 {t)dt] (7) 
<=i Jo 

and Eq. © can be recast as 

%{t) :=%(t)^-(t). (8) 

The introduction of time-dependent "aging" rates presents one possibility 
to describe the non-Markovian effects. It is not unique though. A different 
scheme follows by defining [112]: 

ipij(r) := Ptj(r)if)j(T) (9) 

with Y2iPij( T ) = 1- The physical interpretation is as follows: The process 
stays in a state j for a random time interval characterised by the probability 
density "0j(t). At the end of this time interval it jumps into another state 
i with a generally time-dependent conditional probability Py(r). Evidently, 
any process of the considered type can be interpreted in this way. By equating 
Eq. (jBJ) and Eq. (jHJ) and taking into account ipj{r) := —d^j(r)/dT one can 
deduce that the approach in Ref. [110] can be reduced to that in Ref. [112] 
with the time-dependent transition probabilities 



1 This observation which can be traced back to Ref. [Ill] can be rationalised as follows. 
Let us consider a sojourn in the state j characterised by the survival probability $j(r). 
The corresponding residence time interval [0, r] can be arbitrarily divided into two pieces 
[0, ti] and [ti,t]. If no memory effects are present, then $j(r) = $j(r — ri)$j(Ti). 
The only nontrivial solution of this latter functional equation which decays in time reads 
$j(r) = exp(— 7jT), with > 0. 
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and with the non-exponential probability densities ipj(r) which follow as 

^i( r ) = 7i(i")exp[-/ T 7i(*)dt] with 7j (r) := EiM r )- 

The description of non-Markovian effects with the time- dependent tran- 
sition probabilities Pij(r), is rather difficult to deduce immediately from the 
sample trajectories of an experimentally observed random process The 
same holds true for the concept of time-dependent rates. These rates cannot 
be measured directly from the set of stochastic sample trajectories. On the 
contrary, the RTD ipj(r) and the time-independent p^j (with pa := 0) can 
routinely be deduced from sample trajectories measured, say, in a single- 
molecular experiment [113]. Fig. El renders these assertions more obvious. 
The study of the statistics of the residence time- intervals allows one to obtain 
the corresponding probability densities ipj(r) and, hence, the survival prob- 
abilities $j(t). Furthermore, the statistics of the transitions from one state 
into all other states allows one to derive the corresponding conditional prob- 
abilities pij. From this primary information a complementary interpretation 
of experimental data in terms of time-dependent rates fcy(r) can readily be 
given as 

h 3 {r) = - Pij V d * )l , (11) 

if one prefers to use this language to describe the non-Markovian effects. 
Moreover, the description with a constant set p^ provides a consistent ap- 
proach to construct the stationary realizations of £(t), and thus to find the 
corresponding averaged quantum evolution [114]. 



2.3 Averaging the quantum propagator 

The task of performing the noise- averaging of the quantum dynamics in Eq. 
(J2J) can be solved exactly due the piecewise constant character of the noise 
[34,49]. Let us consider the time-interval [to,t] and take a frozen realization 
of £(t) assuming k switching events within this time- interval at the time- 
instants ti, 

t < h < t 2 < ... < t k < t. (12) 

Correspondingly, the noise takes on the values ^j ,Cju ■■■j^jk * n ^ ne time se- 
quel. Then, the propagator S(t,t ) reads 

S(t,t ) = e~^^fc^* - * fc )e -i£ ^ fc - 1 ^* fc ~* fc_1 \..e -iC ^<'^* 1 ~* ^ . (13) 

Let us assume further that the process £(t) has been prepared in the state 
jo at to- Then, the corresponding k— times probability density for such noise 
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realization is 



p k(t,j k ,t k ;t, jk _ 1} t k -i; •••;Cii)*i|Oo>*o) - ®j k (t-tk)*Pj k jk-i( t k-tk-i)---*Pj 1 j (ti-t ) 

(14) 

for k 7^ and P (^j ,t ) = $j (t — t ) for = 0. In order to obtain 
the noise-averaged propagator (S(t\t ,j )) conditioned on such nonstation- 
ary initial noise preparation in the state jo one has to average (fT3|) with 
the probability measure in (j!4)l (for k = 0,oo). This task can be easily 
done formally by use of the Laplace-transform [denoted in the following as 
A(s) := J °° exp(— sr)A{r)dT for any time-dependent quantity A(r)]. The 
result for (S(s\t , j )) = exp(—sr)(S(t + r\t , j ))dr reads [110,114] 

(S(sKj )) = J2(Ms)[I-B(s)}- 1 ) , (15) 

i 

where the matrix operators A(s) and B(s) reads in components 

POO 

Ms) := 6 U / Mr)e- {s+iCmr dr, (16) 
J o 

and 

/■oo 

B kl (s):= / ^ kl (r)e-^ c ^ T dr , (17) 
Jo 

correspondingly, and / is the unity matrix. 

To obtain the stationary noise averaging it necessary to average (|15jl in 
addition over the stationary initial probabilities pf Q . The averaging over the 
initial distribution alone is, however, not sufficient to arrive at the stationary 
noise-averaging in the case of non-Markovian processes since the noise re- 
alizations constructed in the way just described still remain non-stationary. 
This principal problem is rooted in the following observation. By prepar- 
ing the quantum system at to = in a nonequilibrium state p(0), the noise 
will be picked up at random in some initial state £j with the probability 
pf Q (stationary noise). However, every time when we repeat the preparation 
of the quantum system in its initial state, the noise will already occupy a 
(random) state £,- for some unknown, random time interval r* Q (setting a 
clock at to — sets the initial time for the quantum system, but not for the 
noise, which is assumed to start in the infinite past, cf. Fig. |2l where £j = £i 
at to — 0)- Therefore, in a stationary setting a proper averaging over this 
unknown time r* is necessary. The corresponding procedure implies that the 
mean residence time (r,) is finite, (r,) ^ oo, and yields a different residence 
time distribution for the initial noise state, ipjfr); namely, it is obtained 
as i/)j^(r) = < &j(t)/(tj) [115]. Only for Markovian processes where $j(t) 
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is strictly exponential, does ipf\r) coincides with ipj{r). Using this ipf\r) 
instead of ipj{r) for the first sojourn in the corresponding state and for the 
time-independent set Pij, the noise realizations become stationary [114-116]. 
The corresponding expression for the quantum propagator averaged over such 
stationary noise realizations has been obtained in Ref. [114], cf. Eqs. (25), 
(29) therein. In a slightly modified form it reads 



(S(s)) = (S(s)) statlc - (C(s) - A(s)[I - PD(s)}~ l PA^ 



fa) ' 



where (S(s)) sta tic is the Laplace-transform of the statically averaged Liouville 
propagator 

k 

Pj f = lim^oo Pj (t) are the stationary probabilities which are determined by 
a system of linear algebraic equations [114, 116], 

tr? R "tr (20) 

and P is the matrix of transition probabilities p^ ( "scattering matrix" of the 
random process £(£)). Furthermore, the auxiliary matrix operators C(s) and 
D(s) in (j!8)l read in components: 



Cki{s):=5 kl e"( s+i£ ^ ^(T')dr'dT (21) 
Jo Jo 

and 

Dki{s) := 5 kl / Mr)e- {s+lCmr dr . (22) 



This very same averaging procedure can be applied to any system of linear 
stochastic differential equations. 



2.3.1 Kubo oscillator 



A prominent application of this general procedure is the noise-averaging pro- 
cedure for the Kubo phase oscillator [48,108]; reading 

X(t) = ie[ffl]X(t) ■ (23) 

This particular equation emerges in the theory of optical line shapes, in the 
nuclear magnetic resonance [32,48], and for single molecule spectroscopy [60]. 
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It appears also naturally within our approach, see below, where X(t) corre- 
sponds to a diagonal matrix element of the evolution operator of a quantum 
system with fluctuating eigenenergies. In the context of the stochastic the- 
ory of spectral line shapes [32,48,60], e[£(£)] in Eq. (|2~3|) corresponds to a 
stochastically modulated frequency of quantum transitions between the lev- 
els of a "two-state atom", or transitions between the eigenstates of a spin 
1/2 system. 

The spectral line shape is determined via the corresponding stochastically 
averaged propagator of the Kubo oscillator as [48] 

J(w) = - lim Re[S(-iu + rj)} . (24) 

7T r?— >+0 



Note that the limit rj — > +0 in Eq. (j24j) is necessary for the regularisation 
of the corresponding integral in the quasi-static limit (rj) — > oo. Upon 
identifying with —e k in Eq. (JTBJ) we end up with 

<*(•)> = E^-E^^K (25) 
^ s-ie k ^ {s- te k ) 2 (r k ) 



l-ipi(s-iei) / 1 \ 1 - j) n (s - ie n ) p s n 

^ s — iei \I - PD(s) ) im mn s-ie n (r„) ' 



nl,m 



where D nm (s) = 5 nm ip m (s — ie m ). 2 The corresponding line shape follows 
immediately from Eq. (|23|) by virtue of Eq. (|24|). This result presents a 
non-Markovian generalisation of the pioneering result by Kubo [48] for ar- 
bitrary iV-state discrete Markovian processes. The generalisation consists 
in allowing for arbitrary non-exponential RTDs ipk{r), or, equivalently, in 
accordance with Eq. (fTT|) also for time-dependent transition rates kij(r). 
This generalisation was put forward originally in Ref. [114] for a particular 
case, pf = (Tj)/ $2fc( T fc)j which corresponds to an ergodic process with uni- 
form mixing (meaning that in a long time run each state j is visited equally 
often) . 

Let us next apply this result to the case of two-state non-Markovian noise 
with pi2 = P21 = 1 and pf 2 — { T i,2) /[( r i) + ( r 2)]- Then, Eq. (J25J) yields after 

2 Note that the formal solution of another prominent problem of the first order relax- 
ation kinetics with a fluctuating rate, pit) = — r[£(t)]p(i) follows immediately from (|25[) 
upon substitution ej — > iTj, see in Ref. [105] for some nontrivial non-Markovian examples 
and the corresponding discussion. 
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some simplifications: 

(S( s )\ = y 1 ^ | ( £ i ~ e ?) 2 

x t 1 ~Ms -ie 2 )} , 2g . 

1 - V>i(s - «ei)^ 2 (s - ze 2 ) 

With ()26j) in one obtains the result for the corresponding spectral line 
shape which is equivalent to one presented recently in Ref. [117] by use of 
a different method. It is reproduced within our treatment as a particular 
two-state limiting case. Moreover, in the simplest case of Markovian two- 
state fluctuations with ^1,2(3) = 1/(1+ ( T i,2) s ) and with zero mean, (£(£)) = 
(Ti)ei + (r 2 )e 2 = 0, this result simplifies further to read 



In (|ZZ|), <r= v/(£ 2 (t)) = |e 2 -ei| v /(ri)(r 2 )/((r 1 ) + (r 2 )) denotes the root mean 
squared (rms) amplitude of fluctuations. Moreover, \ = v '/2 + icr sinh(6/2) 
is a complex frequency parameter, where u = 1/ (ti) + 1/ (t 2 ) is the inverse 
of the autocorrelation time of the considered process 3 which possesses the 
autocorrelation function (£(£)£(£')) = cr 2 exp(— u\t — t'\). Furthermore, b = 
ln((ri)/(r 2 )) = In |e 2 /ei| is an asymmetry parameter. The spectral line shape 
corresponding to (J2"7)) has been first obtained by Kubo [48]. It reads [48,64], 



a 2 u 



Moreover, the expression (|2"7j) can be readily inverted into the time domain. 
It is crucial that the corresponding averaged propagator (S(t)) of Kubo os- 
cillator [70], i.e., 

(S(t)) = e-*\coa(yfi*^t) + -f== sin( y/^~^t)] , (29) 
L V o- - X J 

is complex when the process is asymmetric, 6 / 0. This correlates 
with the asymmetry of the corresponding spectral line shape, I{—oS) ^ 
Derived in a different form [118] (for a two-state Markovian process with a 
nonzero mean and in quite different notation) an expression equivalent to 
pnjl is used in the theory of single-molecule spectroscopy [118-120]. For 
a symmetric dichotomous process (with 6 = 0) Eq. (|2HJ) reduces to the 
expression (6.10) (with u = 0) in Ref. [47]. 

3 Note that throughout this work v is the inverse of the autocorrelation time. It is equal 
to the sum of two rates. 
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2.3.2 Averaged dynamics of two-level quantum systems 
exposed to two-state stochastic fields 

The outlined non-Markovian stochastic theory of quantum relaxation can 
be exemplified for the instructive and relevant case of a two-level quantum 
system, reading 

H(t) = Ek\l)(l\ + E 2 \2}(2\ + ^(f)(|l><2| + |2)(1|), (30) 

which is driven by a two-state non-Markovian stochastic field = ±A 
with identical RTDs, i/ji(t) = ^ij) = ip(r). This stochastic field causes 
(dipole) transitions between two states, |1) and |2), and is zero on average. 

This archetype model exhibits a very rich behaviour. In particular, it 
allows one to study the problem of quantum decoherence of a two-state atom 
under the influence of two-state "i/y a " noises exhibiting long range time- 
correlations with a power law decay (for ip(r) possessing a long-time algebraic 
tail, t/>(t) oc l/r 3_a , < a < 1) [121,122]. It thus presents a prominent 
problem of general interest. Moreover, it relates to activities for solid state 
quantum computing [123]. It is convenient to express the Hamiltonian (j3*Uj) 
in terms of Pauli matrices, a z := |1)(1| - |2)(2|, a x := |1)(2| + |2)(1|, a y := 
i(|2)(l| - |1)(2|) and the unity matrix J, 

H(t) = he a z + ^Ht{t)* x + ^{E, + E 2 )I, (31) 

where eo = {E\ — E 2 )/h. Then, the dynamics of the density matrix of the 
quantum two-state quantum system can be given as pit) = | [I+J2 i=x y z &i{t)(Ti] 
in terms of a classical spin dynamics (with components a.i{t) = Tr(p(i)<Xj)) 
in a magnetic field. This latter dynamics evolves on a Bloch sphere of unit 
radius (i.e., the (scaled) magnetic moment is conserved, 4 \(r(t)\ = 1). It 
reads, 

& x (t) = -€ (Ty{t), 

(T y {t) = e a x (t) - £(t)<r z (t), (32) 

The above theory can readily be applied to a noise averaging of 3-dimensional 
system of linear differential equations ()32|) over arbitrary stationary real- 
izations of £(£)• After some algebra, the following result is obtained [124] 

4 This means that each and every stochastic trajectory runs on the Bloch sphere. The 
averaged Bloch vector (a(t)) becomes, however, contracted | ((?(£)) I < 1, because (ai(t)) 2 < 
(af(t)). Thus, the averaged density matrix (p(t)) is always positive in the considered 
model, cf. [47] , independent of the particular model used for the stochastic driving £ (t) . 
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for the Laplace-transformed averaged difference of populations (a z (t)) = 
(pn(t)) — (p22(t)) with the initial condition cr 2 (0) = l,a x>y (0) = , i.e., 
the state "1" is populated initially with the probability one: 

where 

A zz (s) = el[l-i,(s)]{(tt 2 -s 2 )(l-$(s + m)$(s-in)) 

- 2itts[i)(s + iO) -^(s -ifl)]} 

- A 2 s 2 [l + ^(s)][l-^(s + zfi)][l-^(s-^)], (34) 
5 2Z (s) = ^[l-^(s)][l + ^(s + ifi)][l + ^(s-ifi)] 

+ A 2 [l +^(s)](l -i){s + iVL)i){s-iti)), 

and := a/cq + A 2 . Furthermore, r is the mean residence time between the 
field's alternations. Note that for the considered initial condition, (a x (t)) = 
(<J y (t)) = for all times. For e = the result in (J33 )) -(|53 J) reduces to one 
for Kubo oscillator (|26|) with identical ^1,2 (t). Moreover, for the Markovian 
case, if)(s) = 1/(1 + ts), Eq. reduces to 

/- ,u s 2 + 2z/g + z/ 2 + e 2 

S 3 + 2us 2 + (A 2 + £ 2 + ^ + A2l/ ' W 

where v = 2/r is the inverse autocorrelation time. This latter result repro- 
duces the result for the averaged populations (pu(s)) = (1/s + (a z (s)))/2 
and (p 22 (s)) = (1/s — (cr z (s)))/2 in [50,51]. The same result can 
also be reduced from a more general solution for the Markovian case with 
an asymmetric field of non-zero mean [54]. It possesses several remark- 
able features. First, the asymptotic difference between populations is zero, 
(a z (oo)) = lim s _ > o(s(cr z (s))) = 0. In other words, the steady state popula- 
tions of both energy levels equal 1/2, independently of the energy difference 
Hsq. One can interpret this result in terms of a "temperature" T a of the 
(quasi-) spin system. This spin-temperature is formally introduced by us- 
ing for the asymptotic distribution an Ansatz of the Boltzmann-Gibbs form, 
(pnn(oo)) = exp[-E n /k B T a )/Y, n ex-V[-En/kBTo]. Then, 5 

T a := P ^ (36) 



5 This is a standard definition of the temperature of a spin subsystem in nuclear mag- 
netic resonance and similar areas [59] . It is used also to introduce the parlance of formally 
negative temperatures. 
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for two-level systems. In accord with this definition, the result of equal 
asymptotic populations, (p 2 2(°o)) = (pn(oo)) = 1/2 can be interpreted in 
terms of an infinite temperature T a = oo. This constitutes a general find- 
ing: a purely stochastic bath corresponds to an apparent infinite temper- 
ature [40,41]. For this reason, such stochastic approaches to describe the 
relaxation process in open quantum systems is suitable only for sufficiently 
high temperatures k B T ^> h\e \ [40,41]. An asymmetry of unbiased stochas- 
tic perturbations does not change this conclusion, see in Ref. [54] . Moreover, 
the relaxation to the steady state can be either coherent, or incoherent, de- 
pending on the noise strength and the autocorrelation time. In particular, 
an approximate analytical expression for the rate k of incoherent relaxation, 
(Pu{t)) = [1 +exp(— kt)]/2, has been obtained in a limit of small Kubo num- 
bers, K := A/u <C 1, which corresponds to a weakly coloured noise [47,49]. 
This analytical result reads [50,51,54] 

k = (37) 

and exhibits a resonance feature versus v at v = 6q. A similar such resonance 
feature occurs also in the theory of nuclear magnetic resonance for a weakly 
coloured Gaussian noise [59]. Note that in [54] this notable result has been 
obtained for asymmetric fluctuations of the tunnelling coupling possessing a 
non-vanishing mean value ^ 0. This corresponds to a quantum particle 

transfer between two sites of localisation which are separated by a fluctuating 
tunnelling barrier. A related problem with the inclusion of dissipation has 
been elaborated in [67] within a stochastically driven spin-boson model. 

Yet another interesting solution can be obtained for (a x (s)) with the 
initial condition reading a x (0) = 1. The Laplace transform of the solution is 
obtained as 

s 2 + A 2 2A 2 effi 2 A xx (s) 

~ s(s 2 + Q2) TS 2 {S 2 + Q2) 2 £ > (38) 



where 



A xx {s) = [i-4>(s)][i-i>(s + in)}[i-4>(s-in)}, (39) 

B xx {s) = e 2 [l + ^(s)][l - 4>(s + - 4>(s - itl)} 

+ A 2 [l -ip(s)}(l - 4>(s + in)4>(s -in)). 



The physical relevance of this solution (}38|) is as follows. In a rotated quasi- 
spin basis, i.e., a x — ► oy, a z — > ov, a y — > oy, the considered problem 
becomes mathematically equivalent to the problem of the delocalisation of a 



19 



quantum particle in a symmetric dimer with the tunnelling coupling €q under 
the influence of a dichotomously fluctuating energy bias £(£). Therefore, it 
describes the corresponding delocalisation dynamics and, in particular, allows 
one to determine whether this dynamics is coherent or incoherent, depending 
on the noise features. 

For the Markovian case Eq. (JHHj) reduces to 6 



Note that the denominators in Eq. (|35jl and Eq. |40|) are different. 7 In a 
more general case of asymmetric Markovian noise, the corresponding denom- 
inator is a polynomial of 6th order in s, see in [54]. In the considered case 
of symmetric noise it factorises into the product of two polynomials of 3d 
order, those in the denominators of Eq. (|35|) and Eq. (|40|). Thus, for a gen- 
eral initial condition the relaxation of a two-level quantum system exposed 
to a two-state Markovian field involves generally 6th-exponential terms. As 
a matter of fact, this seemingly simple, exactly solvable model can exhibit an 
unexpectedly complex behaviour even in the simplest Markovian case of a 
coloured noise driving. However, for certain initial conditions, as exemplified 
above, the general solution being a fraction of two polynomials of s simplifies 
to the results in Eq. (J35J) and Eq. (jjpj) . 

In a general case of non-Markovian noise, the analytical solutions in Eqs. 
(|3*3*|) and (|3*5|) can be inverted to the time domain numerically by use of a 
numerical Laplace inversion procedure such as the one detailed in Ref. [125]. 

2.4 Projection operator method: a primer 

Next we shall introduce the readers, following Ref. [126], into the projection 
operator technique. We elucidate this scheme by addressing an example that 
is of physical interest in its own right. 

6 The corresponding dynamics also exhibits a resonance feature versus v in a certain 
limit [56]. 

7 A remarkable feature is, however, that both corresponding secular cubic equations 
have the same discriminant, -D(A, v 1 eq) = 0, separating the domains of complex and real 
roots. Hence, the transition from a coherent relaxation (complex roots are present) to an 
incoherent relaxation (real roots only) occurs at the same values of noise parameters, inde- 
pendently of the initial conditions. The corresponding phase diagram separating regimes 
of coherent and incoherent relaxation (judging from the above criterion) has been found 
in [56]. It must be kept in mind, however, that the weights of the corresponding exponen- 
tials are also of importance for the character of relaxation process. These weights naturally 
do depend on the initial conditions. 




(40) 
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Let us consider the somewhat more general dynamics, 



da(t) 
dt 






and let us pose the question: How can we extract a single, closed equation 
for the evolution of cr z (t) (without any approximation) for an arbitrary time- 
dependence of the parameters governing the driven quantum dynamics? The 
use of a projection operator method provides an elegant way to solve this 
problem [126]. The key idea is to project the whole dynamics onto a corre- 
sponding subspace of reduced dimensionality by using a projection operator 
V with the idempotent property V 2 = V. In the present case, the choice of 
this projection operator follows naturally as: 



(42) 



The use of this projection operator allows one to split the whole dynamics 
into the "relevant" one c?o(t), and a remaining, "irrelevant", //(£), part, re- 
spectively; i.e., a(t) = Va(t) + (1 — V)a(t) := a (t) + jl{t) by applying V 
and the complementary projection operator 1 — V to Eq. (|41jl . From the 
resulting system of two coupled linear equations for <7o(i) and jl(t), i.e., 





^ = VB(t)a (t)+VB(t)fI(t), 

^ = (l-V)B{t)ff (t) + (l-V)B(t)m, ( 4 3) 
a single integro-differential equation for aa(t) follows, reading: 
da (t) 



(Jf = VB(t)a (t) 

+ VB{t)Texp (J dr(l - V)B{r) \ (1 - V)B(t')a (t')dt' 

+ VB{t)Texp Qf dr{\ - V)B{t) \ jl(0) . (44) 

The exponential matrix operations in ^4*]) can be done explicitly without any 
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approximation, yielding the exact closed equation for a z (t) [126], reading 8 



is introduced which is a functional of the time- varying parameter e(i). The 
projection of the entire dynamics onto the some subspace typically entails 
memory effects. Put differently, a nonlocality in time emerges for the reduced 
space dynamics. Moreover, an explicit dependence on the initial conditions 
in the "irrelevant" subspace is necessarily present. 

3 Two-state quantum dynamics in periodic 
fields 

Let us illustrate the practical usefulness of the exact equation (|4*5j) by its 
application to a quantum dynamics occurring in strong time-periodic fields. 
Towards this goal, we consider a quantum two-state tunnelling system, where 
the two states |1) and |2) correspond to the two sites of charge localisation 
(i.e. we work in the "tunnelling" representation) and A = const corresponds 
to the tunnelling matrix element. This charge dynamics is driven by a pe- 
riodic electric field of frequency Q which results in a periodic modulation of 
the energy bias between two localised states e(t) = eo + Acos(flt) [3]. We 
assume that the particle is prepared initially on the site "1" at t — 0, i.e. 
cr z (0) = 1 and a x (0) = o~ y (0) = 0. Further, one assumes that the frequency 
of external field is rather high, Q ^> A and we consider the correspondingly 
averaged dynamics (cr z (t))n using the high-frequency decoupling approxima- 



tion (cos[((t,t')]cr z (t'))n « (cos[C(t, t')]) n (a z (t')) Q [83]. Using Eq. (gED and 



the well-known identity exp(iz sin 6) = Yl^L-oo Jn{z) exp(m6>) (J n (z) denotes 
the Bessel function of the first kind) the high-frequency approximation in Eq. 

8 Within the path-integral approach, the same equation can be derived from a non- 
interacting blip approximation (NIBA) result of the dissipative spin-boson model [132] by 
putting formally therein the strength of the system-bath coupling to zero. Astonishingly 
enough, the NIBA turns out to provide the exact result for this singular limit of zero- 
dissipation. 




+ A(t) sin[C(t, 0)K.(0) + A(t) cos[C(t, 0)K(0) . 



(45) 



In Eq. (JISJ), the time-dependent phase 




(46) 
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(1351) yields 



ft 

2 



(<7,(t)) n = —A / T(t - t')(a z (t')) n dt' (47) 
J o 

with the kernel r(t) = Xl^L-oo ^(-^A^) cos [( e o + ?if2)t]. This latter equation 

fit 



can be solved by the use of the Laplace-transform. For a z (s) := f Q exp(— st) 
(crz(t))ndt one obtains 

&z(s) — - r>, a TFt • (48) 

1 ^ ^ Z^n=-oo s 2 + ( eo _ n Q)2 

From this relation now follow some key-results. 

3.1 Coherent destruction of tunnelling (CDT) 

The formal inversion of the result in Eq. ()48jl into the time domain reads 

(v z (t))n= cjexpiiujjt), (49) 

j=0,±l,±2,... 

where Uj = isj > are the poles Sj of Eq. pHj) . From the quasi-periodic 
character of the driven dynamics it follows that all these poles lie on the 
imaginary axis in complex conjugated pairs. Therefore, uj_,j = —Uj and 
c_j = c*. Although there appears an infinite number of poles, only few of 
them contribute significantly in the regimes of interest. 

Let us consider the case of symmetric two-level system (TLS), eo = 0. 
Then, the approximate solution reads (only the term with n = in the sum 
in Eq. (|48|) contributes significantly in the high-frequency limit A/Q — > 0) 

(vz{t))n = cos(At un t), (50) 

where A tun = A\J (A/Q)\ is the renormalised tunnelling frequency. When 
the amplitude A of the high-frequency driving is chosen to obey J (A/fl) =0, 
the tunnelling dynamics is brought (within this high frequency approxima- 
tion) to a complete standstill. This constitutes the celebrated phenomenon 
of coherent destruction of tunnelling [127, 128] which attracted much atten- 
tion and generated many applications over recent years, see e.g. in [3,4] and 
references therein. 
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3.2 Driving-induced tunnelling oscillations (DITO) 



Let us consider now the case of a large energy bias eo 3> A. In the ab- 
sence of driving, the particle remains essentially localised on the site "1", 
<j z (t) ~ 1, as can be deduced from the well-known exact solution a z {t) = 
[eg + A 2 cos(y / eQ + AH)]/[el + A 2 ]. When, however, a high-frequency driving 
ft 3> A is applied such that the resonance condition nft = a/cq + A 2 « e 
is approximately fulfilled, large amplitude tunnelling oscillations in Eq. (|50|1 
can be induced with a tunnelling frequency A tun = A\J n (A/Q)\. This phe- 
nomenon of driving-induced tunnelling oscillations (DITO) - being opposite 
to CDT has been revealed in Refs. [87,129,130]. It has recently been verified 
and observed experimentally [130]. This DITO-phenomenon is illustrated in 
Fig. E] for an "exotic" 5-photon (n = 5) resonance case where the precise 
numerical solution of driven TLS dynamics using Eq. (J4~T|) and the approxi- 
mation in Eq. (}50|) with A tun = A\J n (A/Q)\ are plotted for the following set 
of parameters: A = 0.1, e = 20, A = 24, ft = 3.9 or ft = 4.0. For ft = 3.9 
the dynamics is almost localised exhibiting small-amplitude oscillation - cf. 
dotted line near a z (t) = 1 which is barely visible in Fig. El a but becomes 
clearly seen in Fig. El b due to a better resolution on a different scale. A 
relatively small change of the periodic field frequency chosen to match the 
resonance condition 5ft = eo induces large-amplitude tunnelling oscillations 
which are nicely described on the long-time scale by the approximation in Eq. 
(joT?|) with A tun — A\J 5 (A/fl)\. It cannot be distinguished from the precise 
numerical solution in Fig. El a. The frequency of these oscillations is con- 
trolled by both the bare tunnelling frequency A and the field amplitude A. It 
is worth mentioning that DITO seems to be close in spirit to the famous Rabi 
oscillations [131] (an interpretation given in Ref. [130]), but are in fact by 
no means identical with those. Rabi oscillations correspond usually to a par- 
ticular case of a small-amplitude, A <C ft, resonant driving, ft = a/A 2 + eg. 
For €q = 0, the corresponding problem is equivalent (in a rotated quasi-spin 
basis, b x — > ay, a z — > oy, b y — > a y >, and with different initial conditions) 
to the resonant dipole excitation of a two-state atom with eigen-frequency 
Ld = A. In such a case, the frequency of Rabi oscillations ojr is determined 
approximately by the driving amplitude, i.e. ujr « A [3,161]. This presents 
the most remarkable, characteristic feature of Rabi oscillations. The DITO- 
frequency presents rather a driving- renormalised tunnelling frequency as in 
the case of CDT. The coarse-grained character of the result in Eq. (|5T)|) is 
illustrated in Fig. El b on a short-time scale in comparison with the precise 
numerical solution of the driven dynamics. This latter one exhibits step-like 
transitions with a number of oscillations on each step. The number of oscil- 
lations corresponds to the number of emitted (absorbed) photons. With the 
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Figure 3: Driving- induced tunnelling oscillations. 

increase of n, the "steps" become longer and sharper. In order to make a 
further "step" in the transfer of population the two level system awaits for 
the next portion of n photons to be emitted, or absorbed to match the res- 
onance condition nQ = eo (a quasiclassical interpretation of the numerically 
observed step feature). 

4 Dissipative quantum dynamics in strong 
time-dependent fields 

4.1 General formalism 

Without loss of generality we consider a iV-level, driven quantum system 
characterised by a time-dependent Hamilton operator Hg(t) and which inter- 
acts Vsb with a thermal bath characterised by a Hamilton operator H B - The 
system-bath interaction is assumed here to be generally also time-dependent. 
It is characterised by the Hamilton operator Vs B (t) which depends both on 
the (relevant) variables of the system of interest and on the thermal bath 
variables. The total Hamiltonian H(t) thus reads 

H(t)=H s (t)+V SB (t)+H B . (51) 

The dynamics of the density operator pit) of the total system is then gov- 
erned by the corresponding Liouville-von-Neumann equation, cf. Eq. (|T|). 
Furthermore, the reduced density operator of interest is obtained by per- 
forming a partial trace of p(t) over the bath variables, i.e. ps{t) = Trgp(t). 
The average < A > of any operator A which depends on the variables of 
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the system of interest can be calculated as the corresponding trace over the 
system variables, i.e., < A >= Trg(pg(t)A). The reduced density operator 
ps{t) , which also depends on the initial preparation scheme, thus contains 
all the necessary information required to describe the time-evolution of the 
system of interest. The main task consists in obtaining a tractable closed 
equation of motion for ps(t). This can be achieved by applying to pit) a 
properly chosen projection operator II, which projects the whole dynamics 
onto the subspace of the considered quantum system, thereby accounting 
indirectly for the "irrelevant" bath variables, i.e. ps{t) = Itp(t). A proper 
choice for the projection operator with the idempotent property, II 2 = II, is 
II := p B Tr B [22-24], where ps = Gxp(—PH B )/Z B is the equilibrium density 
operator of the bath; Z B = Tr B exp(— (3H B ) is the corresponding partition 
sum, and (3 = l/{k B T) denotes the inverse temperature. Then, pit) can 
identically be split as p(t) = p B ® Psif) + r ?( i 0> where r](t) = Qp(t) represents 
a cross-correlation term. Here, Q := 1 — II is the complementary projection 
operator with the properties Qli = HQ = 0, Q 2 = Q. By applying II and Q 
to the Liouville- von- Neumann equation for p(t), two coupled linear operator 
equations for ps{t) and r] (t) can be obtained, respectively, which in turn yield 
a single closed equation for psit) after having eliminated the part r](t). The 
formally exact equation for the reduced density operator, thus reads 

Ps (t) = -iL s (t) Ps (t) - [ T(t, t')p s it')dt' + h[t), (52) 

Jo 

where 

T(M') = TT B [L SB it)S s+ Bit,t')QL SB it')p B } (53) 
denotes the memory kernel. In Eq. (15*31 . 

S s+B (t,t') = Texp{-2 / [L s {t) + L b + QL SB ir)]dT} (54) 

Jt' 

is a Liouvillian propagator. Furthermore, Ls(t)(-) = [Hs(t),(-)]/h, L B {-) = 
[H B , i-)]/h, LsBit)i~) = \YsB(t), i-)]/h are the corresponding Liouville opera- 
tors, where H$it) := Hsit)+ < VsB(t) >b is the renormalised Hamiltonian 
of the dynamical system and Vss(t) := Vs B (t)— < Vs B it) > B is the corre- 
spondingly re-defined system-bath coupling. 9 Moreover, J (t) in Eq. (p)2*j) 

/ (t) = -zTr B (L SB (t)S s+B (t, 0)/i(0)) (55) 

9 This issue deserves to be commented on in further detail: The generalised quantum 
thermal forces acting on the system from the bath should be on average unbiased. This 
implies that the thermal average < ... >b'= Trs(ps...) of a properly defined system-bath 
coupling, < VsBit) >b-= Trs(/9BVsB(t)), should be zero, i.e. < Vss(i) >b= 0. For this 
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constitutes the initial correlation term, sometimes also termed "initial value 
term" . 

Note that the GME (f52|) - (f55|) is still exact in the subspace of the quan- 
tum system for a quantum evolution started at to = 0, i.e. no approxima- 
tions have been invoked so far [26,27,133]. Generally, a reduced quantum 
evolution contains some dependence on the initial conditions /i(0) in the ir- 
relevant subspace. We note, however, that for a factorised (uncorrelated) 
initial preparation p(0) = ps <S> Ps(0) (yu(0) = 0) this initial correlation term 
vanishes identically, i.e. Io(t) = 0. This standard class of initial preparations 
will be assumed in the following. 

4.1.1 Weak-coupling approximation 

In the second order approximation with respect to the system-bath coupling 
VsB(t) (the so termed weak-coupling limit) one sets LsB(t) — > in Ss+B(t, t'), 
Eq. (JHH). Moreover, let us assume a factorising form of the system-bath cou- 
pling Vss(t) = \ J2 a K a{t)laia + h.c. where 7 Q denote some system opera- 
tors, £ a are the bath operators, and K a (t) are the coupling strength functions. 
The complete set j a is assumed to be closed under commutation relations: 
[7a, 7/?] = J2s e a/3S A /5 with e af 3s being some structural constants defining a 
corresponding Lie algebra with generators The Hamiltonian Hs(t) is 
represented as a linear superposition 



in the corresponding algebra. For N-level quantum systems the following set 
of operators is conveniently used. It is given (here a := (n,m)) by the set of 
operators 7„ m := \n)(m\ (they correspond to the elements of corresponding 
Liouville space). Here the the ket- vectors \n) provide an orthonormal vector 
basis with the scalar product (n\m) = 5 nm in the corresponding Hilbert space 
of the considered N-level quantum system. The representation of the system 
Hamiltonian in this discrete basis reads 



with H nm (t) = H^it). It is evident that any quantum system with a discrete 
number of states can be represented in this way. The system-bath coupling 

reason, the systematic, mean-field like contribution < VsB(t) >s of the thermal "force" 
should be separated from the very beginning and be included in Hs(t) without change of 
the Hamiltonian of the total system. Obviously, this can always be achieved. This formal 
renormalisation is always assumed in the following (with "tilde" omitted when applicable). 




a 




(56) 



ran 
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can be chosen in the form 

Vss{t) = 

nm 

with K mn (t) = K* nm {t) and £ mn = £} nm . Moreover, the dissipative operator 
kernel in Eq. (p)2*|) reads in the given approximation: 

r(^,^)(') = ^ 1 K nn' {t ) K mm > (t ) ^ K nn ' mm i (t t ) 7nn') S (t , t )j mm i (■) 
n,n' ,m,m' 

~ Kn'nm'mit ~ t') Inn' , S(t,lf)(-)^ mm ' }, (58) 

where 

S{t,t') = Texp{-i [ L s {r)dT} (59) 
Jt' 

is the Liouville evolution operator of the physical system under considera- 
tion. Note that it does include the external, time-dependent field influences 
exactly. Moreover, 

K nn i mm '(t) := — < inn'{t)imm< >B= K m > mn > n( — *) > (60) 

is the autocorrelation tensor of the thermal force operators £ n n'(t) '■= e %HBi l h 
inn'£~ %HBt l h - An expression formally similar to Eq. (J58j) has been obtained 
first, for a particular case of a spin 1/2 system (and for a time- independent 
system-bath coupling) in Ref. [24]. For the reduced density matrix, p nm (t) := 
(n\ps(t)\m) the following generalised master equation (GME) follows: 



L nmn> m >{t)p n , m/ {t) - ^ / {t,t')p n , m ,{t')dt', (61) 



where L nmn > m >(t) = \[H nn >(t)S mm > - H m , m (t)5 nn >] is the Liouville superopera- 
tor, written in the supermatrix representation and the corresponding memory 
kernels read 

r nmn 'm' it , t') = | Knki{t)Kkn'(t')K n k'kn'(t ~ t')Uk'k(t,t')U* nm ,(t,t')(62) 
kk' 

~l~ ^k'm{t)K m 'k(t ^■^mk'km 1 

(t-t')U nnl (t,t')U* k , k (t,t') 

K nk >(t)K m 'k(t')K k , nkm ,(t - t')U k , n r(t,t')U^ k (t,t') 
- Kk'm{t)K>kn'{t)Kk'mkn'{t ~ t) U nk (t, t ) Uy m , (t, t ) 

28 



where U mrn /(t,t') := (m\T exp{ — | f* H s (T)dT}\m') is the evolution operator 
of the considered quantum system in the Hilbert space. This result presents 
the most general form of weak-coupling GME in arbitrary external fields. 
Generalised master equations of a similar form have been derived before, by 
making use of different methods and in different notations in Refs. [64, 66, 
82] 10 . The kernel (I62|) satisfies two important properties which must be 
strictly obeyed; these are (i): 

^nmn'm'it-, — ^mnm'n'i^i ^ ) 

(imposed by the requirement that ps{t) must be Hermitian, ps(t) = pjj(£)), 
and (ii): 

^ ^ ^ rinn'm' ifi t ) 
n 

(conservation of probability, Tr $ ps{t) = 1 for all times). 

This driven GME (|6ip presents our "working horse" that will be used 
frequently for the discussion of various applications detailed below. 



4.1.2 Markovian approximation: 

Generalised Redfield Equations 

The integro-differential equations ([61 )1 - ([62 )1 are, notably, non-local in time, 
i.e., they describe a so-termed "non-Markovian" quantum dynamics. This 
non-locality in time makes their practical use rather cumbersome. A cor- 
responding Markovian approximation, which renders a description that is 
local in time, is therefore of great use in practice, if it can be justified on 
physical grounds. There are several ways to obtain such a Markovian ap- 
proximation. The most popular one is to perform a back propagation, i.e., 
p s (f) = S ,_1 (*,i / )Ps(*) + 0(k 2 ), in the kernel of GME making use of the 
Liouville evolution operator S(t,t') of the dynamical subsystem. The corre- 
sponding master equation for the reduced density matrix, which constitutes 
the generalisation of the well-known Redfield equations [135] to the case of 
driven, open quantum systems reads: 

Pnm(t) — ~i L nmn r m '(t)p n r m '(t) — Rnmn'm' {t) Pn'm' (t) j (63) 

n'm' n'm' 

10 For a periodic driving, the field influence on the relaxation kernel can alternatively 
be taken into account applying the corresponding Floquet basis for periodically driven 
quantum dynamics. This has been done in Refs. [84,95]. Our approach is valid, however, 
for arbitrary time-dependence. 
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with the explicit time- dependent relaxation tensor 



K n l(t)Kkk<(t')K n ikk>(t — t)Ulk(t,t)U*, k ,(t,t)8 mm i 
Klm(t)Kkk'(t')KZ uklk (t - t')U mlk {t,t')Ui k ,{t,t')8 nn i\ 

K nn >(t)K W (t')K*, nk , k (t - t')U m/k (t,t')U^ kl (t,t') 

Km'm(t)K kk >(t')K m > mkk i(t — t')U nk (t,t')U*, kl (t,t')\dt' . 

(64) 

This relaxation tensor satisfies two important relations, namely, R nm n'm' (t) — 
Kanm'n'if) (Ps(t) is Hermitian), and J2 n R nnn'm'(t) = (imposed by the 
conservation of probability). Notable, the upper limit of integral in (|6"4"j) is 
here the actual evolution time t (instead of oo) - this feature softens already 
the well-known problem which relates to a possible violation of positivity at 
initial time scales of the quantum evolution [136] for Redfield equations for 
certain initial conditions 11 . Upon neglecting (setting to zero) the influence of 
external time-dependent fields in the relaxation tensor, by using the basis of 
eigen-states of time-independent H s , and setting t — ► oo in (jMJ), we recover 
the commonly known form of the Redfield relaxation tensor. 

It must be stressed that the physical nature of the thermal bath operators 
was up to now still not specified. Those can be either be of bosonic, fermionic 
nature, or also describe a spin bath [134] . The corresponding autocorrelation 
tensor (|fiO"|l has to be calculated for every particular microscopic model. We 
next address within this methodology several physical applications. 

11 This problem can be resolved by the so-called slippage of the initial conditions, see 
in [137-139]. Moreover, within the weak-coupling approximation the effect of dissipation 
should be consistently taken into account to the second order of the system-bath coupling 
only; i.e., in the solutions of Redfield equations (rather than in the relaxation kernels only). 
The dissipation-induced frequency shifts (i.e. the Lamb shifts at T = 0) should also be 
very small (compared to the corresponding eigen-frequencies of quantum evolution in the 
absence of dissipation). Otherwise, the theory needs to be renormalised. Notwithstanding 
these essential restrictions, the Redfield equations provide one of the most widely used 
tool to describe open quantum systems in many areas of physics and physical chemistry 
[3,4,9,59,140-143]. 



Rnmn'm' (0 ^ 
kk' 



kk' 

t 
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5 Application I: Quantum relaxation in 
driven, dissipative two-level systems 

Let us consider a two-level quantum system with time- dependent eigenenergy 
levels 12 , 

H s (t) = [E{ 0) + E 1 (t)]\l){l\ + [E^+E 2 (t)}\2)(2\, (65) 

which is coupled to a bath of independent harmonic oscillators, possessing 
the spectrum {ui\}, i.e., 

H B = J2^x(b{h + l), (66) 

A 

where b\ and b\ are the bosonic creation and annihilation operators, respec- 
tively. The interaction with the thermal bath causes the relaxation transi- 
tions between the eigenstates of the dynamical system ( "longitudinal" inter- 
action). Such transitions are absent otherwise and thus require either the 
emission, or the absorption of bath phonons. The interaction is chosen to be 
of a bi-linear form, reading 

VfcB = |(|l>(2| + |2)<l|), (67) 

with the influence of the bath being presented by a random force operator, 
£12 = £21 = £, 

I = + (68) 

A 

Here and elsewhere below, the (real-valued) coupling constants k's are in- 
cluded into the fluctuating force £. From a phenomenological perspective, 

12 These levels can correspond, e.g., to spatially separated localisation sites of a trans- 
ferring (excess) electron in a protein [11]. If such electronic states possess very different 
dipole moments (the difference can reach 50 D [144]), an external time-dependent electric 
field will modulate the energy difference in time due to the Stark effect. Such an electric 
field dependence of the electronic energy levels can be very strong [144, 145]. A large 
modulation of the local electric field can be induced, e.g., due attachment/detachment of 
an ATP molecule/products of its hydrolysis. A substantial shift of the electronic energy 
levels can then be induced [146]. In a simple setting, the corresponding modulation of an 
energy level can be modelled by a two-state Markovian process [147]. The chemical source 
of driving force can also be substituted by a direct application of a stochastic electric 
field [147,148]. This latter possibility has been demonstrated experimentally for some ion 
pumps [148]. 
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the considered model represents an analogue of the model in Sec. 12.3.21 where 
a classical random force assuming two values is replaced by a quantum oper- 
ator force which possesses a Gaussian statistics. Moreover, a possible time- 
dependence of the energy levels is assumed here. 

The corresponding correlation function of this Gaussian quantum stochas- 
tic force K(t) := K 122 i(t) =< ^12(^)^12(0) > is complex- valued. It reads 
explicitly 

K{t) = — [ J(cu)[coth(-^-)cos(iut)-ism(ujt)]duj. (69) 
27r J 2k B T 

with the bath spectral density given by 

J ^) := §rE*-^). (70) 

A 

Upon extending to negative frequencies u < 0, it is convenient to formally 
define J(— uj) := — J{ui). 

The complex nature of this bath correlation function is crucial for the 
establishment of thermal equilibrium at the finite temperatures. The ap- 
plication of GME (JBTj) . (HI2J) to the present case yields a closed system of 
generalised master equations for the level populations populations p n (t) '■= 
Pnn(t),n = 1,2, reading, 

p\(t) = - [ {w 12 {t,t') Pl {t')-w 2l {t,t')p 2 {t')]dt', 
Jo 

p 2 {t) = [ [w 12 (t,t') Pl (t , )-W2i(t,t , )p 2 {t'))dt', (71) 
Jo 

with memory kernels 13 

w 12 (t,t') = 2Re[K(t-t')exp(ie (t-t')+i((t,t')), 

w 21 {t,t') = 2Re[K(t-t')exp(-ie Q (t-t')-i((t,t% (72) 

where e = (£{ 0) - Ef ] )/h and £(t, t') is a functional of time-dependent driv- 
ing, Eq. (gHJ), with e(t) = [E^t) - E 2 (t)]/h. We shall assume that e(t) fluc- 
tuates (either randomly, or periodically in time) around a zero mean value. 
In order to obtain the quantum relaxation averaged over the fluctuations of 
e(t) one needs to perform a corresponding stochastic averaging of the GME 
([71)1 . For an arbitrary stochastic process e(i), this task cannot be carried 
out exactly any longer; consequently one must resort to some approximation 
scheme(s). 

13 One can immediately see that if K{t) would be real- valued, then the forward and 
backward rate kernels were always equal. This would mimic the situation of an infinite 
temperature as it was elucidated in Sect. 2 for a classical stochastic bath. 
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5.1 Decoupling approximation for fast fluctuating en- 
ergy levels 

If the characteristic time scale for e(t) fluctuations r e is very small in compar- 
ison with the characteristic system relaxation time scale r r , i.e. r e <C r n then 
one can use a decoupling approximation by averaging for (pi !2 (t)) e ; namely, 
(exp(±i((t, t')pi^{t')) e « (exp(±i((t,t')) e (p^ 2 (t')) e . For fast fluctuations of 
the energy levels the relaxation dynamics then follows (pi, 2 (t)) t with fast, su- 
perimposed small-amplitude fluctuations whose amplitude diminishes when 
the ratio r e /r r <C 1 becomes smaller. A subsequent Markovian approxima- 
tion for the averaged dynamics yields a master equation description of the 
form: 

(pi(t)) e = -(W 12 (eaMPi(t))e + {W 21 (e )) e {p 2 (t)) e , 

(P2(t))e = (W 12 (e ))e(pi(t)) e -(W 21 (e )) e (p 2 (t)) t (73) 

with the averaged transition rates reading [64], 

/oo 
e^n(ou)J(u])I(e -uj)duj, 
■oo 

/oo 
n(u)J(uj)I(e -uj)du . (74) 
-oo 

Here n(u) = l/[exp(hw/ (ksT)) — 1] is the Bose function, and I(uj) is the spec- 
tral line shape of a Kubo oscillator X(t) = ie(t)X(t) (see in Sect. 2). From 
Eqs. (JZHJ), (j71j) one can immediately deduce that in the absence of fluctua- 
tions, where I (to) = 5(u>), the thermal equilibrium, pi(oo) = e^^^^p^oo), 
is attained independently of the specific model for J(u>) for arbitrary tempera- 
tures T. Moreover, the thermal detailed balance condition, p 2 (oo)W 2 i(e ) = 
Pi(oo)Wi2(eo) is obeyed always with the thermal bath temperature T. In 
other words, the temperature T a of the considered TLS, defined through Eq. 
()36|) coincides with the temperature of the thermal bath T, T a = T. This is 
in a sharp contrast to the stochastic bath modelling in Sec. 12.3.21 where we 
found that pi(oo) = ^2(00) and T a = 00. Furthermore, one can see that the 
thermal equilibrium at the bath temperature T becomes violated either by 
periodic, or by stochastic nonequilibrium fluctuations. This also implies that, 
generally, T a 7^ T. Put differently, either periodic, or stochastic (thermally 
nonequilibrium) fields drive the system out of the thermal equilibrium with 
the thermal bath. This fact lies at the heart for the emergence of a diver- 
sity of interesting and often counter-intuitive nonequilibrium effects which 
we shall address next. 
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Figure 4: (a) The averaged relaxation rate of TLS, ro(eo), in units of 
nl/(h 2 Qo), is depicted versus the noise amplitude a (in units of f2 ) f° r 
the averaged energy bias e = 0.4f2 and the thermal bath temperature 
T = 0.25htt /k B and 7 = 0.051V (b) The effective temperature of TLS, 
T a , in units of M^o/Zcg versus the noise amplitude a (in units of Qq) for the 
same set of parameters of TLS and the environmental temperature T. At 
a ~ 0.6f2 , where the TLS is maximally cooled, the lower level is populated 
with the probability being close to one. On the contrary, for a 1.4O the 
upper level becomes populated with the probability being close to one; i.e., an 
almost complete inversion of populations occurs. The model assumptions are 
well justified for the coupling constant kq <C 0.05M)q such that r (eo) <C 7. 



5.1.1 Control of quantum rates 

A first application is the manipulation of the transition rates by many or- 
ders of magnitude by use of a rapidly fluctuating, discrete stochastic fields 
[62,64-66]. This scheme becomes feasible when the spectral density of the 
bath J{uj) is sharply peaked around some vibrational frequencies. The effect 
can be demonstrated for a quantum Brownian oscillator model of the bath: 
It corresponds to a single quantum vibrational mode VLq which acquires a 
frictional broadening width 7 due to a bi-linear coupling to other environ- 
mental vibrational modes 14 . The corresponding spectral density assumes the 

14 A fast (on the time scale r r of system relaxation) equilibration of this single mode 
with other vibrational modes is assumed. This imposes an important restriction t^ 1 := 
(Wi2(eo)) e + (W / 2i(eo))e 7 which can always be justified by a proper tuning of the 
coupling constant k . Furthermore, the broadening of vibrational spectral lines in molec- 
ular systems 7 exceeds typically 7 > 5 cm -1 (in spectroscopic units) which corresponds 
7 > 10 12 Hz in units of the frequency. The considered relaxation transitions must conse- 
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form [5]: 



J H = tt /. .o 02^ , ro ' ( 75 ) 



/1 2 (u 2 - n 2 ) 2 + 4 7 2 w 2 

Let us start out by considering first a control scenario of quantum relaxation 
by use of a symmetric dichotomous Markovian noise e(t) = ±cr with I (us) 
given by (J2EJ), where ei )2 = ±0". For the case 15 v <C cr, this spectral line shape 
consists of two sharply shaped peaks, located at uj = ±cr and possessing the 
width v. For v <C 7, which is typically the case, this latter broadening can 
be neglected. Then, I(uj) « |[5(o; — cr) + 5(c<j + cr)], and the averaged rates 
simplify to read 

(W 12 (e )) e « ^ 12 (e + cr) + ^ 12 (e -cr)] 

4/tg7^ / (eg - cr)e fc a T n(e - cr) 
ft 2 v[(eo-cr) 2 -fi 2 ] 2 + 4 7 2 (eo-cr) 2 

(e + cr)e fc a T n(e + cr) \ 

[( eo + (T )2_ fi 2 ]2 + 472(eo + (T)2y ) I J 

(W 2 i(e )) e « i[^ 21 (e + cr) + ^ 21 (eo-cr)] 

4^7^ / (eo - o")n(e - cr) 



+ 



ft2 V[( eo -cr) 2 -fi 2 ] 2 + 4 7 2 (eo-cr) 2 
(eo + cr)n(e + cr) 



[(e + cr) 2 -fi 2 ] 2 + 4 7 2 (eo + cr) 



From Eq. (|7fij) it follows that if the quantum transition frequency eo matches 
the vibrational frequency of the medium Qq an increase of energy fluctuation 
size cr (induced by local electric field fluctuations in the medium) from zero 
to some finite value cr ^> 7 can drastically reduce the relaxation rate r (e ) = 
(Wi 2 (eo)) e + (W / 2i(eo))e! even a practical blockade of relaxation transitions 
can occur [62,64-66]. On the contrary, in the case of a frequency mismatch 
between eo and Qq one in turn can dramatically enhance the rate of relaxation 
transitions by tuning the noise amplitude cr appropriately [66] , see in Fig. HI 
a. 



qucntly occur more slowly, e.g., an electron transfer can occur on a msec time scale [11]. 

15 This case presents a realistic situation for experimental realizations of molecular 
systems since a significant stochastic perturbation with an energy exceeding one fc^T, 
ha ~ 25 meV (at room temperatures), corresponds in the units of frequency a ~ 4 x 10 13 
1/s. The frequency v of large amplitude bistable conformational fluctuations of molecular 
groups is typically much smaller. 
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5.1.2 Stochastic cooling and inversion of level populations 

A second effect relates to the blockage of the rate for backward transitions 
(W / 2i(eo))e relative to the forward rate (W / i2(eo))e- This can cause a stochastic 
cooling of the TLS, where the temperature T a becomes smaller than the 
temperature of the environment, i.e., T a < T. This interesting phenomenon 
is demonstrated with Fig. |H(b). Similar in spirit, although different in the 
physical mechanism is the laser cooling (of the nuclear degrees of freedom) 
as it has been studied both, theoretically and experimentally for polyatomic 
molecules [149]. 

Moreover, for a > Q a noise-induced inversion of steady state averaged 
populations takes place, i.e., for a sufficiently small positive energy bias e the 
higher energy level becomes more populated. This constitutes the third very 
important effect under discussion, see Fig. IH b, where this pumping effect 
is interpreted in terms of a negative temperature T a . In other words, the 
considered nonequilibrium noise of a sufficiently large amplitude is capable 
to pump quantum particles from the lower energy level to the higher one. 
This provides a possible archetype for quantum molecular pumps driven by 
nonequilibrium noise. 

This inversion of population can be accompanied by cooling. Namely, 
the ensemble of TLSs becomes first effectively cooled and only then heated 
up (through, formally, T a — oo to T a — — oo) until the inversion of popu- 
lation occurs - cf. in Fig. 01 b. For this pumping mechanism to work, an 
inverted transport regime [184] is necessary; i.e., a regime where the static, 
unfluctuating forward rate becomes smaller with the increasing energy bias 
after reaching a maximum at e max . In the present model, this maximum 
is located in the neighbourhood of fl . More precisely, e max corresponds 
to the maximum in the difference between the forward and the backward 
rates, rather than to the maximum of the forward rate alone. The inver- 
sion happens for a > e max and a sufficiently small energy bias e . A similar 
mechanism has been proposed in Ref. [68] within a spin-boson modelling of 
electron tunnelling in proteins, see below in Sec. El driven by nonequilibrium 
conformational fluctuations, e.g., utilising energy of ATP hydrolysis. 

The underlying mechanism seems quite general. Indeed, the inversion of 
populations occurs whenever the difference of averaged rates (AW(eo)) e '■= 
(Wi2(eo)) e — (W / 2i(eo)) e becomes negative, (AW(eo)) e < 0, for a positive bias 
e > 0. In the discussed limiting case, (AW(e )) e ~ \{AW(t + a) + AW(e - 
a)] with AW(— e) = — AW(e). Therefore, when a exceeds e max , where 
AW(e) achieves a maximum and 4-AW(e) < for o > e max , the averaged 
difference of forward and backward rates becomes negative (AW(eo)) t < 0, 
for a positive energy bias eo > 0, i.e. an inversion of populations takes place. 
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In application to the quantum transport in a spatially extended system, a 
similar effect results in the noise-induced absolute negative mobility [69], see 
below in Sec. The existence of the static current-voltage characteris- 
tics with a negative differential conductivity part is important for the latter 
phenomenon to occur. 

5.1.3 Emergence of an effective energy bias 

The fourth important effect, the onset of which can be seen already in the dis- 
cussed archetype model, is rooted in a possible asymmetry of the unbiased on 
average fluctuations. Namely, let us consider the symmetric quantum system, 
eo = 0, with asymmetric dichotomous fluctuations of the energy levels with 
zero mean, see in Sec. 12.3.11 Since in this case, the averaged propagator of 
the corresponding Kubo-oscillator is complex-valued, Im (exp[i({t,t')]) t ^ 0, 
it can be readily seen from Eq. (|72|) after invoking the decoupling approx- 
imation that {wi2(t,t')) e — (w2i(t,t')) t 7^ even if eo = 0. This means 
that an effective asymmetry emerges. Moreover, the above difference is pro- 
portional also to Im Kit — t') ^ 0. If the autocorrelation function of the 
thermal bath, K(t), were real (like for a stochastic bath), then no asym- 
metry between the forward and backward rates could emerge in principle. 
Therefore, the discussed asymmetry does emerge due to a subtle interplay 
of the equilibrium quantum fluctuations of the thermal bath and nonequilib- 
rium classical fluctuations of the energy levels, both of which are unbiased 
on average. Here is rooted the origin of quantum dissipative rectifiers put 
forward in Ref. [70,102]. The very same effect can also be deduced from 
Eq. (HH), since the corresponding spectral line I(uj), cf. Eq. (}2"%j) . is asym- 
metric. Yet, ultimate insight is achieved in the slow-modulation limit of the 
Kubo oscillator, K e :— ajv ^> 1, like in Eq. (j7(ij) where the mean forward 
and backward rates are the static rates Wi 2 (e) and W 2 i(e) averaged over the 
energy bias distribution, pit\ 2 ) = (ti2)/(( t i) + ( r 2)), correspondingly, i.e., 
(W 12 (0))e = E W f(«ei) and (W 21 (0)) e = E^pfe^rfe). 

For several applications of quantum transport in spatially extended sys- 
tems our driving-induced breaking of symmetry leads to a rectification cur- 
rent in tight-binding Brownian rectifiers [70,102], see in Sec. 

5.2 Quantum relaxation in strong periodic fields 

The considered strong nonequilibrium effects are present as well in the case 
of a fast periodic driving, e(t) = A cos(flt + ip ) , with a static phase ipo which 
is uniformly distributed between and 2tc. Then, the corresponding spectral 
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line shape form I(oj) is = Y^L-oq Jn(A/ty8(u} — nQ), where J n (z) is 
the Bessel function of the first kind. 

The rate expressions (|76jl take on the form 

00 A M ^] 

(W 12 (e ))e = Yl fi^e^^nfa-ntyJfa-nto), 

n=— 00 

00 . 

(W 21 (e )) e = E J 2 n(^)n(e -nQ)J(e -nQ). (77) 

n=— 00 

Such an expansion of the transition rates as a sum over different emission 
(absorption) channels with n emitted (absorbed) photons with the corre- 
sponding probabilities p n = J%(A/fl) is similar to one used by Tien and 
Gordon in a different context [3,150]. For the averaged relaxation rate the 
above expression yields the same result as in [82] where the principal possibil- 
ity to regulate quantum relaxation processes in condensed molecular systems 
by strong periodic external fields has been indicated. Moreover, the inver- 
sion of populations by periodic driving takes also place for the above model 
J(oo) and some properly adjusted parameters of the periodic driving. For 
a periodically driven spin-boson model (see below in Sec. 6) and a strong 
system-bath coupling, this latter effect has been theoretically predicted and 
described in Refs. [83,89] (see also [3] for a review and further references). 
For the case of a weak system-bath coupling, the inversion of populations in 
the spin-boson model has been shown in Ref. [87]. 

5.3 Approximation of time-dependent rates 

If the external field varies sufficiently slow on the characteristic time-scale, r^, 
describing the decay of the kernels in Eq. (fTTj) . an adiabatic approximation 
of time-dependent rates that follow the time-variation of the energy levels 
can be invoked; i.e., 

Pi(t) = -W 12 (e(t)) Pl (t) + W 21 (e(t))p 2 (t), 

p 2 (t) = W l2 {e{t)) Px {t)-W 2X {e{t)) V2 {t), (78) 

where W 2 i(e) = n(e)J(e) and Wi 2 (e) = exp(he/k B T)W 2 i(e) are the static 
rates. For discrete state noise, this approximation holds whenever (tj) ^> r^. 
Then, the corresponding rates describe a discrete state stochastic process 
and the averaging method of Sec. El can be applied. In this limiting case, the 
corresponding Laplace-transformed averaged populations can be given in ex- 
act analytical form. In the considered case this corresponds to the averaging 
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of a Kubo oscillator with an imaginary frequency. The corresponding aver- 
aged solution can be analytically inverted into the time domain in the case of 
two-state Markovian fluctuations. This solution is generally bi-exponential. 
Two limiting cases can be distinguished which can be classified by the Kubo 
number K\y of the rate fluctuations, i.e. by the product of the variance of 
the rate fluctuations multiplied with the corresponding autocorrelation time. 
In a slow modulation limit (in terms of the rate fluctuations), Ky/ 3> 1 and 
the (ensemble) averaged relaxation is approximately described by a quasi- 
static averaging of the time-dependent solutions with the static, "frozen" 
energy bias randomly distributed. It assumes a bi-exponential form, but can 
be multi-exponential and anomalously slow in a more general case of multi- 
state fluctuations. The opposite limit of fast modulation (in terms of the rate 
fluctuations, Ky/ <C 1) corresponds to the averaged rate description which is 
detailed above and which invokes the decoupling approximation and in addi- 
tion can possibly involve a slow modulation limit in terms of the energy level 
fluctuations, K e ^> 1. The resulting averaged relaxation process remains 
approximately single-exponential. 

In view of the presence of many different time scales the underlying 
physics is nontrivial. Therefore, it is useful to be able to resort to a case 
study where the stochastic averaging can be performed exactly. Such an 
archetypal investigation has been put forward in Ref. [100] and has been 
applied to the stochastic spin-boson model in Refs. [68,101]. 

5.4 Exact averaging for dichotomous Markovian fluc- 
tuations 

By use of the conservation of probabilities, the system of integro-differential 
equations (J7T|) can be reduced to a single equation for the difference of pop- 
ulations a z {t) = pi(t) —p2(t). It reads: 



*z(t) = - [ f(t,t')a z (t')dt' - [ g{t,t')dt' 



(79) 



Jo Jo 



with the integral kernels 



/(*,*') 

g(t,t') 



fo(t,t') cos[e (t 
g (t,t') sm[e (t 



+ C(M')L 




where 



/o(M') 

9o(t,t') 



f (t-t') = ARe[K(t-t% 
g (t-f) = -A lm[K{t-t')}. 




39 



The kernel f(t,t') in Eq. (J75J) denotes a stochastic functional of the driv- 
ing field e(t) on the time interval [t',t] (with times later than t') whereas 
a z (t') is a functional of the dichotomous Markovian process (DMP) for the 
times prior to i! . The task of stochastic averaging of the product of such 
functionals, (f(t,t')a z {t')), can become very difficult, see [151]. However, 
in the case e(t) = aa(t), where a(t) = ±1 is symmetric DMP with unit 
variance and autocorrelation time r c = 1/V this task can be solved exactly 
by referring to a theorem by Bourret, Frisch and Pouquet [152] (for a dif- 
ferent proof of this remarkable exact result, see also Ref. [153]). It reads: 
(f(t, t' + T)a(t'+r)a(t')a z (t')) = (f(t, t' + r)){a{t ! + r)a{t')){a z {t ! )) + (f(t, t'+ 
r)a(t' + r))(a{t')a z (t')) for r > 0. By passing to the limit r — ► and us- 
ing the characteristic property of the DMP, namely that a 2 (t) = 1 (without 
averaging), the above relation yields an important corollary [100]: 

(f(t,t')a z (t')) = (f(t,t'))(a z (t')) + (f(t,t')a(t'))(a(t')a z (t>)) . (82) 

This result is beyond the decoupling approximation, given by the first term 
in the sum. The result for the cross-correlation function (a(t)a z (t)) is more 
intricate. The equation of motion for this cross-correlation function can be 
obtained due to a theorem by Shapiro-Loginov [51, 154], reading, 

j t (a(t)a z (t)} = -u(a(t)a,(t)) + • (83) 

Use of this relation in turn generates an integro-differential equation for 
{a(t)a z (t)), where the problem of decoupling of (a(t)f(t,t')a z {t')) emerges. 
It can be solved in the same way as in Eq. (|82|). namely 

(a(t)f(t,t')a z (t')) = (a(t)f(t,t'))(a z (t')) 

+ (a(t)f(t,t')a(t'))(a(t')a z (t')) . (84) 

All these averaged functionals, like (f(t,t')), (a(t)f(t,t')), (f(t,t')a(t')), 
(a(t)f(t,t')a(t')) can be expressed in terms of the averaged propagator of 
the corresponding Kubo oscillator S^(t — t') = (expficr L a(r)dr]) given 
in Eq. (j2Uj) with x = i//2 (zero asymmetry), and its derivatives S^(t) := 
ij^°)(f) [100,101]. Applying the general results in Eqs. (jH2), Eq. jgg) 
and (JH1|) to Eq. (|T§j) yields a closed system of two integro-differential equa- 
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tions [100,101]: 



j t Mt))= - / {S^(t-t')f (t-t')cos[e (t-t')){a z (t')} 

- SW{t - t')f (t - sin[e (t - t')}(a(t')a z (t')) 

+ S^(t-t')g (t-t')sin[e (t-t')]\dt', (85) 



j t (a(t)a s {t)) = - v{a(t)a z (t)) 

t 



S&(t - 0/o(* - cos[6 (t - f)]<«(*>*(0> 
+ - t')f (t - sin[e (t - t')](a z (t')) 

+ S {1 \t - t')g {t - t') cos[e (t - t')))dt' . 



A subsequent Markovian approximation for Eq. (|53|) then yields [100]: 

j t {° z {t)) = -T ((T z (t)) -T 1 (a(t)a z (t)) - r , 
-(a(t)a z (t)) = -Y 1 {a z (t))~{u + Y 2 ){a{t)a z {t))-r 1 (86) 



with 



coth ( — — — ) J(uj)Ik(e — lu)cIlu 



2k B T / 

PCX 

r k = J(uj)I k (e - uj)duj , (87) 



oo 



where h(uj) = {—oj/a) k I{uj) and J (a;) is given in Eq. (j28|) with ei j2 = 
±a. It can be shown that all known limiting cases are reproduced from this 
remarkable result. 

For the case of weakly coloured noise K t 1 (i.e. the fast modulation 
limit of the Kubo oscillator), the spectral line I(oS) becomes a Lorentzian 
with the width D = a 2 jv. The same result holds true in the white noise 
limit a — > oo, v — > oo, with D = const and K e — > 0. In these limits Ti is 
negligible small, I\ « 0, and the relaxation is described by the averaged rate 
r . Precisely the same result re-emerges also for white Gaussian noise e(t) 
with the noise intensity D. The spectral line becomes narrower when v 
increases, - this constitutes the celebrated motional narrowing limit of NMR 
[32,33] -, and approaches zero when v — > oo (with a kept constant). Such 
infinitely fast fluctuations have no influence on the considered rate process; 



41 



the field-free description is thus reproduced with the thermal equilibrium 
being restored. 

In the slow modulation limit of the Kubo-oscillator (K € 3> 1), ~ 
h{uS) ~ + a) + 5(u — a)] and h{uj) rs \[8{uj + a) — S(u — a)], with 

the corresponding line widths neglected. In this case, the approximation of 
time-dependent fluctuating rates following adiabatically to the energy levels 
fluctuations becomes justified. The relaxation is generally bi-exponential 
with the two rates given by 



Here, T± = coth[h(e ±a)/2kBT]J(e ±cr) is the relaxation rate in the quasi- 
static limit and To = (T + + r_)/2. Furthermore, if v 3> To, we have Ai ~ v 
(note that the corresponding exponent exp(— Ait) contributes, however, with 
a very small weight), and A2 ~ T (with a weight which approximately equals 
one): This in turn implies that the relaxation is practically single exponential 
with rate T and corresponds to the fast modulation limit in terms of the 
fluctuating rates. 

6 Application II: Driven electron transfer 
within a spin-boson description 

Let us proceed with an application of our general theory to the celebrated 
driven spin-boson model [3]. This model is of special importance since it 
describes a large variety of physical phenomena [1,2], such as relevant aspects 
of electron transfer in molecular systems. 

6.1 Curve-crossing problems with dissipation 

The simplest case of a two-state donor acceptor electron transfer reaction 
can be considered as a curve-crossing problem within the description of two 
diabatic electronic states |1) and \2) with electronic energies V\(x) and V^x) 
that depend on a nuclear reaction coordinate x [8,9,157-159] (cf. Fig. EJ). 
Namely, after separating nuclear and electronic degrees of freedom within the 
Born-Oppenheimer approximation, the electron tunnelling process coupled 
to the nuclear dynamics (modelled by the reaction coordinate x) can be 
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Figure 5: Sketch of the diabatic electronic curves. Note that two crossing 
points occur in the presence of different curvatures [160] . 



described by the following Hamiltonian 



P 

2M 



+ V 1 {x,t) 



+ ~£A(t)(|l)<2| + |2><l 



The time-dependent electronic curves in Eq. 



1 



2M 



+ V 2 (x,t) 



|2)(2| 



Vi, 2 (x,t) = -Mni 2 (x±x /2y±h6 /2-d 1>2 £(t), 
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(90) 



can generally possess different curvatures in the parabolic approximation 
with minima energetically separated by he and separated by a distance x 
(the tunnelling distance). Moreover, such electronic states generally possess 
electric dipole moments d\ }2 (their coordinate dependence is neglected) and 
thus the discussed energy levels will generally become dependent either on 
the stochastic microscopic fields of the environment, or on an externally 
applied electric field S(t). The corresponding time-dependence can likewise 
reflect also some nonequilibrium conformational dynamics. Moreover, the 
coupling or tunnelling matrix element A(£) can also parametrically depend on 
a nonequilibrium reaction coordinate which generally introduces an explicit 
stochastic time- dependence. The reaction coordinate x is coupled to the 
rest of vibrational degrees of freedom of the environment. This introduces 
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a dissipation mechanism into the tunnelling problem which can be modelled 
by a bilinear coupling of x to the thermal bath of harmonic oscillators [5,20], 

VTii 



H ^= oE — ^-—2* • (91) 



It is worthy to point out that the frequencies f2i )2 of the oscillator x can 
depend on the electronic state. In other words, the relevant vibration modes 
can become either softer, or more rigid depending on the electronic state. In 
the following we neglect this possible effect and assume that Qi = Q 2 — 
but note the studies in Refs. [155, 156, 160] for the more general situation. 
Moreover, one assumes that the reaction coordinate relaxes rapidly (with 
respect to the time scale of electron transfer) into thermal equilibrium with 
the bath of oscillators. Then, a canonical transformation from the "reaction 
coordinate + N bath oscillators to "N+l new bath oscillators" brings the 
original problem into the spin-boson form, i.e., [5] 

11 1 

H{t) = -he(t)a z + -HA(t)cr x + ^x a z ^ c x x\ 

x 

where e(t) = eo — {d\ — d 2 )£(t)/h. The coupling between the quasi-spin 
and boson bath is characterised by a so-called spectral density 16 J(oS) = 
(7r/2) J2i(c\/'™'>'ti\)8{u — &\) [!]• Moreover, we assume for the low frequency 
behaviour an Ohmic-like coupling between the reaction coordinate and the 
environmental vibrational modes (which corresponds in the classical limit to a 
viscous frictional force F = —Qx acting on the reaction coordinate x) which in 
turn yields the effective spectral density J {(J) = (u (^2_o2)2 +4 ^2 7 2 i 7 = C/2M. 
This scheme corresponds to the model of a damped Brownian harmonic os- 
cillator used in Eq. (|7S|) with k q = \J (M7 )A, where A = is the 
reorganisation energy. The coupling strength ( can be related to the dimen- 
sionless (Kondo) parameter a = ^ = | ^ The use of the representation 
of bosonic operators in Eq. ()6.1|) then yields 

Hit) = ^He(t)a z + -HA(t)a x + -a zV (t) J2 K *( b l + 6 >) 

A 

+ J2hw x (b{b x + ±) (93) 



16 This definition is related to the one given in Eq. I|70() by: J(u>) — 2xqJ(uj)/H with 
k\ = x cxy/H/{2rhx6jx), x x = \fh/ \2rhxuJx){b\ + b\). 
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(the "tilde" over uj\ is omitted here). To address formally the most general 
case we assume in addition that the system-bath coupling can be modu- 
lated in time as well, i.e., n\ — > K\r)(t) with some prescribed time-dependent 
function fj(t). 

6.2 Weak system-bath coupling 

Let us consider first the case of a weak system-bath coupling. The corre- 
sponding generalised master equations are obtained by applying Eqs. (JEH), 
(in the representation of 7 nm ) to the considered spin-boson model. This 
yields after some cumbersome calculations by using the quasi-spin basis the 
following GMEs: 

a x (t) = -e(t)<r y {t)- [ T xx {t,t')a x {t')dt'- [ T xy {t,t')a y {t')dt' - A x {t), 

Jo Jo 

& y (t) = e{t)a x {t) - A{t)a z {t) - [ T yx (t,t')a x (t')dt' 



o 



t 

T yy (t,t')a y (t')dt' -A y {t), (94) 



2/2/ V ' " J W V\ U J" 1 " 1 *2/ 

& z (t) = A(t)a y (t), 

with the kernels reading 

r xx (t,f) = V (t) V (t')Re[K(t - t')}Re[U 2 u (t,f) + U 2 2 (t,t% 

r w (t,0 = 77(t)7Kt')Re[X(t - t')}Re[U^(t,t') - U 2 2 (t,t% 

T xy (t,t') = V {t) V (t')Re[K{t - t')}lm[U 2 u (t,t') - U 2 12 (t,t% 

r yx (t,t') = - r)(t)r)(t')Re[K{t - t')}lm[U 2 u (t,t') + U 2 l2 (t,t% (95) 

and the inhomogeneous terms given by 

A x (t) = ll 7 7 (t)7/(t , )Im[ J ft:(*-t / )]Im[U 11 (t,tOU 12 (t,t / )]c?t / , 
Jo 

A y {t) = 2/ r7(^)r/(^')Im[i ; ^(t-^ / )]Re[Ull(t,t / )Ul2(t,t , )]rft , . (96) 
Jo 

The evolution operator of the driven TLS in the absence of coupling, which 
defines the Hamiltonian Hjj(t) of the driven, nondissipative dynamics, is 
denoted by 



U nm (t,t') = (n|Texp 



% -^E D {r)dr 



\m). 
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This propagator enters the above memory kernels; it can be found nu- 
merically from the solution of the corresponding Schrodinger equation for an 
arbitrary time-dependence. Moreover, in the case of a periodic driving, an ex- 
pansion into Floquet modes is conveniently applied, see in [161], and further 
references therein. Other methods, e.g., the use of a Magnus expansion [162] 
are also possible. Due to the unitary quantum evolution in the absence of 
dissipation we have ^(MO = ^11 (MO an d #21 (MO = — ^(MO witn 
det[U nm (t, t')] = 1 for arbitrary time-dependence of e(i). 

Time-nonlocality of the generalised master equations in Eq. (|9*3J) makes 
them difficult to study from a numerical viewpoint. To work with a memo- 
ryless, Markovian description presents, therefore, a pivotal advantage. If the 
dissipation is very weak, this description suffices to capture the main influ- 
ences of dissipation on the driven quantum dynamics, i.e., the emergence of 
an exponential relaxation (and decoherence) described by some small rate 
constants and corresponding, dissipation-induced frequency shifts, i.e. the 
Lamb shifts occurring even at T = 0. Both the relaxation rates and the fre- 
quency shifts are proportional, in the lowest order, to k\. Applying Eqs. (|615|). 
(HID) to the considered dynamics yields the following driven Bloch-Redfield 
equations: 

& x (t) = -e(t)<r v (t) - R xx (t)a x (t) - R xz (t)a z (t) - A x (t), 

& y (t) = e{t)a x {t) - A(t)a z (t) - R yy (t)a y (t) - R yz (t)a z (t) - A y (t), (97) 

& z (t) = A(t)a y (t) 

with the relaxation matrix elements reading 

R xx (t) = Ryy(t) = [ ??(t)?7(t0M^-O][IMM0lMMM0H< 

Jo 

R xz (t) = 2/ r J (t)r } (t')Re[K(t-t')]Re[\] n (t,t')\] l2 (t,t')}dt', 
Jo 

R yz {t) = -2 I ri(t)ri(t')Re[K(t-t , )]lm[U 11 (t,t')U 12 (t,t , )]dt' . (98) 
Jo 

In the common case of a time-independent tunnelling matrix element, i.e. 
A(£) = const and a time-independent system-bath coupling, i.e. rj(t) = 1 
(what is assumed in the following), this result reduces to the driven Bloch- 
Redfield equations derived in Ref. [87]. Note the different signs of A and 
e used throughout this work and in Ref. [87], as well as some other cited 
references. 
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Figure 6: (a) Numerical comparison of the driven Bloch-Redfield equations 
of Ref. [87] (dotted line) and the path-integral GME of Ref. [94] (full line) 
for an oscillatory high-frequency driving Q ^> ujq (data taken from Ref. [87]). 
Both depicted numerical solutions practically coincide within line width. The 
dashed-dotted line depicts a quasi-analytical solution (for details see in [87]) 
of the driven path-integral GME. It captures well the main features of the 
driven dynamics, lacking only some finer details. Time and frequencies are 
measured in units of A -1 and A, correspondingly. The used parameter sets 
are depicted in the figure, (b) Corresponding asymptotic long-time dynam- 
ics: the numerical solution of the driven Bloch-Redfield equations (dotted 
line) is compared with the quasi-analytical solution of driven path-integral 
GME (full line). Both solutions agree well within the width of the small- 
amplitude, driving induced oscillations. The two insets depict the analytical 
results for the rate of averaged relaxation T^(eo) and the difference of asymp- 
totic populations Poo(eo) : = — liim-^ a z (t), respectively. The rate of inco- 
herent relaxation Tr exhibits characteristic resonance peaks, being located at 
multiple integers of the driving frequency f2. These peaks are shifted replicas 
of the dc-driven rate with different weights, i.e. the case with no oscillatory 
forcing acting (i.e. a vanishing driving amplitude s = 0). Thus, a suitable 
chosen static field eo can enhance or suppress the decay of populations. The 
asymptotic population difference P^ exhibits a non-monotonic dependence 
versus the asymmetry eo when combined with a high frequency driving field. 
For appropriate values of bias e , a population inversion takes place (P^ < 
when e > 0, and vice versa). 
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For constant bias e(i) = e , and constant tunnelling coupling A(t) = A, 
U u (t,t') = cos[u {t-t')/2]-i—sm[uJo{t-t > )/2}, 
U 12 (t,t') = -i— sin[uo(t-t')/2] , (99) 

where uq = ^/ + A 2 . Then, the equations reduce to the non-driven Bloch- 
Redfield equations of Ref. [163]. Some different weak coupling master equa- 
tions for the driven spin-boson model have been derived in Ref. [94] using the 
path integral approach. The equation for a z (t) (not shown here) has the form 
of a closed integro-differential equation of rather involved form. In the limit 
of vanishing dissipation it reduces to Eq. (}4"5|) derived within a projection 
operator formalism. 

The numerical equivalence of the our driven Bloch-Redfield equations and 
the weak-coupling integro-differential equation of Ref. [94] has been demon- 
strated in Ref. [87], both by comparison of the numerical solutions of both 
equations for the initial-to-intermediate part of the relaxation time scale 
and by comparison of the numerical solution of the Bloch-Redfield equation 
and an approximate analytical solution of the weak-coupling GME of path- 
integral approach on the whole relaxation time-scale. This numerical compar- 
ison has been performed for periodically driven case, e(t) = — e — s cos(fit), 
for the Ohmic bath with exponential cutoff, J (to) = Anauje^ ul ^ c , where a 
is the dimensionless coupling strength (Kondo parameter) which has to be 
sufficiently small 17 . 

Both approaches agree quite well, cf. in Fig. El The presented approach, 
however, is technically more convenient. The results possess a broad range 
of applications; for example, it allows one to study a mechanism of suppres- 
sion of quantum decoherence by strong periodic fields for a two-level atom 
dynamics in an optical cavity [164]. The investigation of similar mechanisms 
is also of prime importance for the investigation of the quantum decoherence 
in various quantum information processing applications [165]. 

17 An important restriction is: a(A/cjo) 2 hi(^e/^o) <C 1 for lo c 3> ujq. It stems from 
the requirement of the smallncss of the frequency Lamb shift, luq — » Lu r , at T = 0. This 
restriction is most crucial for eo = 0, where o>o = A and uj r = A r rs A[l — aln(w c /A)] ~ 
A exp[— a ln(w c /A)] w A( A)«/(i-«) (f or a < 1, to the linear order in aln((D c /A)). Thus, 
this frequency shift is consistent with the renormalisation in section 6.3. For a large 
asymmetry cq ^> A, the validity range of Bloch-Redfield equations in a becomes broader. 
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6.3 Beyond weak coupling theory: Strong system-bath 
coupling 

Thus far we concentrated on the case of weak-coupling to quantum thermal 
heat bath, or the regime of weak dissipation, respectively. The analytic the- 
ory is, however, not restricted to the case of weak dissipation only. In fact, by 
use of a combination with the method of canonical (unitary) transformations 
one can study the opposite limit of strong dissipation and weak tunnelling. 
To do so, let us consider the spin-boson problem in Eq. (|93J) in the case of a 
strong coupling between the quasi-spin and the bath degrees of freedom. As 
a primary effect, the bath oscillators will become shifted due to this coupling 
to new positions which depend on the spin state. If the tunnelling coupling 
A were absent, then the small polaron unitary transformation [1,11,166-171] 

J7 = exp[i*,J2], R = y £ i £L(Bl-B x ) (100) 

A 

to the new basis of displaced bath oscillators B\ = Wb\U = b\ + ^j-(? z , 

B\ = Wb\U = b\ + ^^cr z and boson-dressed spin states, \n) := W\n) would 
in fact diagonalise the Hamiltonian, solving thereby the problem of finding 
the eigenstates of the total system exactly. For this reason, the corresponding 
canonically transformed basis of phonon-dressed quasi-spin states (polaronic 
states) and displaced bath oscillators is well suited for an approximate treat- 
ment in the case of weak intersite tunnelling and strong system-bath coupling. 
In this new polaronic basis the Hamiltonian reads, 

H{t) = ^e(*)[|l)(l|-|2)(2| 

+ ~ftA(*)( < e k > B |1)(2|+ < e~ k > B |2)(1|) 

+ ^A(t)(V- < e k > fl ]|l>(2| + [e- k - < e~ k >„]|2><1|) 

+ i^nw A (S][B A + l/2)-A//4 (101) 



2 

A 



where 



2Wo w K ' 



is the reorganisation energy. 
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Since < exp[±.R] > B = exp[< R 2 > B /2] = exp [-£)], where D = 
j- J^°[J(lj) coth(j3hw) /uj 2 ]du) and f3 = l/(k B T), the effective tunnelling cou- 
pling, A r := Aexp(— D), between the polaronic states is exponentially sup- 
pressed by the Debye- Waller factor [11,167,168]. For the relevant case of 
Ohmic coupling, J (no) = Amauj exp(— u/u c ), D — > oo and A r — > due to the 
infrared divergence of the corresponding integral. One can attempt to remove 
this divergence by using instead of k,\ in the polaron transformation some 
variational parameters to be determined from the requirement of a minimum 
of the (free) energy of the whole system [170]. An approximate solution of the 
corresponding variational problem by using the Peierls-Bogolyubov-Feynman 
upper bound for the free energy [172,173] leads [170] to a self-consistent equa- 
tion for A r which at T = and for the symmetric case e(t) = reads, 

A r = A exp [ - -!- / . J ( U } .. fa] . (103) 
L 47ri (uj + A r ) 2 J v ' 

Numerically, it can be solved by iterations. An approximate analytical so- 
lution is also available in the limiting case w c > A for a < 1. It yields 

a. 

the celebrated renormalised tunnelling matrix element,i.e., A r = A\J^-j 1 

[1, 2, 170, 174, 175] with uj c = Cu c , where C is some constant which de- 
pends on the precise form of cutoff function in J(oS). In this case, the use 
of the variationally optimised polaron basis allows one to obtain an effective 
Bloch-Redfield description which interpolates well between weak and strong 
dissipation, see for the undriven case the study in Ref. [170]. The correspond- 
ing generalisation of this approach onto the driven case for an intermediate 
coupling strength a < 1 remains yet to be done. Within our approach this 
generalisation is rather straightforward. 

We proceed further with the case of a strong coupling, i.e. a > 1, where 
A r does iterate to zero for any fixed value of u c . This fact indicates the 
famous dissipation-induced localisation transition [1,174,175]. In this case, 
the discussed divergence is not removable; it is real. The polaronic states are 
strictly localised in this case. This is also the feature that causes the local- 
isation phase transition in the dissipative tight-binding model [2,176,177]. 
The third line in Eq. ()101|) presents a (small) time-dependent interaction be- 
tween the dressed system and the bath which can be handled in perturbation 
theory in the lowest order of tunnelling coupling A. Applying the GME (161)1 
to the considered case of an Ohmic bath yields [67, 101] a GME in the form 
of Eqs. (J7HJ), (jSOjl wherein fo(t,t') and g (t,t') assume, however, a distinct 
different form; namely, 

/o(M') = A(t)A(t')exp[-ReQ(t-t')]cos[lmQ(t-t% 

g {t,t') = A(t)A(t r )exp[-ReQ(t-t')]sm[lmQ(t-t% (104) 
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where 



Q(t)= [ dt x f K(t 2 )dt 2 + i\t/h (105) 
Jo Jo 

denotes the doubly-integrated autocorrelation function of the bath, K(t), 
in Eq. 1)69)1 . For A(t) = const the same generalised master equation was 
derived in Ref. [88] using a different approach. It has been derived also in 
Refs. [98] using the path-integral method within the so-called noninteracting 
blip approximation (NIBA). In the case A(t) = const and e(t) = const, it 
reduces to the NIBA master equation of Refs. [178-180]. 

Notably, the driven NIBA master equation is valid for a > 1 at T = 
and eo = (and sufficiently small A <C X/h = 2au c ). It can also be used, 
however, for a < 1 for an asymmetric case, eo ^ 0, and/or for T > 0, where 
the dynamics (in the absence of driving) is incoherent and where A r = 0. 
The parameter domain, where this latter condition is fulfilled, is defined from 
the solution of a (more complicated than Eq. (JlOHjl ) self-consistent equation 
for A r which generally depends on the static bias eo, temperature T, cutoff 
uj c , It can be solved only numerically: in particular, for e ^ and T = 0, 
the renormalised tunnelling coupling vanishes, A r = 0, already for a > 1/2. 
Moreover, even for zero energy bias, eo = 0, the renormalised tunnelling 
coupling vanishes at a sufficiently high temperature, iraksT > HA [170]. 
Even more, for A r ^ 0, the incoherent tunnelling regime holds obviously 
when ksT 3> HA r . Surprisingly, however, for the symmetric situation, eo = 
0, the NIBA master equation turns out to be a very good approximation even 
for arbitrarily small a and T (including coherent dynamics) in the so-called 
scaling limit uj c ^> A with A r fixed. This remarkable fact is rationalised 
within the path-integral approach [2]. Some understanding can be obtained 
by observing that in the limit of vanishing dissipation a — > the NIBA 
master equation is exact and it reduces to the one in Eq. , for the initial 
condition being o~ z (0) = ±1. This, however, amounts to a singular limit 
which must be handled with care. 

6.3.1 Fast fluctuating energy levels 

Let us assume for the following an incoherent quantum dynamics with a 
time-independent tunnelling matrix element A(t) = const. In the case of 
fast stationary fluctuating energy levels the procedure of Sec. El leads (after 
Markovian approximation) to an averaged dynamics in Eq. (|73|) with the 



51 



time-averaged transition rates given by 

(W 12 (e )) e = -A 2 Re / e^- Q « (S(t)) e dt (106) 
* Jo 

1 r°° 

(W 21 (e )) e = -A 2 Re J e '^-Q(t) {S *(t)) e dt, (107) 



where 

(S(t)) £ := (S(t + t ,t ))e = (e*C° eW >, (10? 



is the averaged propagator of the corresponding Kubo oscillator which does 
not depend anymore on the initial time to, or the initial phase of driving. 
These averaged rates can be also given in the equivalent spectral representa- 
tion form, like in Eq. (|74jl. 

poo 

(W 12 (eo))e = -A 2 / FC{u)I{eo-u)du, 

" J -oo 

poo h 

(W 21 (e )) e = -A 2 / e-^FC(co)I(e -co)duj, (109) 

^ J -oo 

where 

1 f°° 

FC(uj) = — exp[iut - Q(t)]dt (110) 

is the Franck-Condon factor 18 which describes spectral line shape due to 
multi-phonon transitions [8,181,182], and I{uS) denotes the spectral line shape 
of the Kubo oscillator, X(t) = ie(t)X(t). The result in Eq. ()106|) is in essence 
the Golden Rule result generalised here to fast fluctuating nonequilibrium 
fields. This fact underlines the generality and importance of the nonequilib- 
rium Golden Rule result which is very useful in many applications. Many 
profound nonequilibrium effects described in this work can be rationalised 
within its framework. The structure of this result has a clear physical inter- 
pretation. Namely, FC(u) in ()110|) is nothing but the spectral line shape of 
a quantum Kubo oscillator with the frequency modulated by the quantum 
Gaussian force £ (t) in Eq. (|fi8j) (in the corresponding Heisenberg representa- 
tion) which has the complex-valued equilibrium autocorrelation function in 
Eq. ([B9*)l. Due to the Gaussian character of the quantum random force, this 
spectral line shape in Eq. ()110j) is expressed merely in terms of the doubly- 
integrated autocorrelation function K(t) and the reorganisation energy term 

18 i.e. the thermally weighted overlap of the wave functions of displaced quantum oscil- 
lators. 
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in Eq. (jlU5j) . Due to the equilibrium character of quantum fluctuations, 
FC(lu) possesses a symmetry property, FC(— to) = e"^FC(w), which is 
enforced by the thermal detailed balance condition. It holds independently 
of the form of the bath spectral density J{uj) [1]. Thus, the thermal equilib- 
rium for localised energy levels 19 , pi(oo) = e~ tu °^ kBT p2(oo), holds always in 
the absence of nonequilibrium fluctuations of the energy levels. Furthermore, 
by splitting £ into a sum of two arbitrary statistically independent compo- 
nents (two subsets of quantum bath oscillators), £ = £1 + £2 one can show 
that FC(u) can exactly be represented as a frequency convolution of the 
corresponding (partial) Franck-Condon factors FC\(u) and FC2(w) [2,8], 
namely, 

/oo 
FC l (oj)FC 2 (e - uj)du. (Ill) 
-00 

Such frequency convolution can be generalised to an arbitrary number of par- 
titions. The nonequilibrium Golden Rule in Eq. (jl()9j) presents an additional 
frequency convolution with the spectral line shape I(uj) of the nonequilib- 
rium Kubo oscillator which corresponds to a generally non-Gaussian and 
nonequilibrium stochastic force. does now no longer possess the above 
symmetry imposed by thermal detailed balance. Thus, the violation of the 
thermal detailed balance condition by the nonequilibrium fluctuations lead 
generally to intriguing nonequilibrium effects described in Sec. and below. 
It is important to notice that the localised states can be stabilised by strong, 
fast oscillating periodic fields [183] and the Golden Rule description is gener- 
ally improved for such fields [97]. This latter fact can be readily understood 
from the representation of the (quantum) stochastic force as a sum of statis- 
tically independent components. Namely, if A r = 0, due the interaction with 
a subset of oscillators, the addition of an interaction with further oscillators 
cannot enhance A, r . It will work always in the direction to make the effective 
tunnelling coupling smaller (when A r 7^ 0), thus improving the perturbation 
theory in A. Replacing equilibrium oscillators with a fast fluctuating field 
does not change this trend. 

6.3.2 Exact averaging over dichotomous fluctuations 
of the energy levels 

An exact averaging of the NIBA master equation of the driven spin-boson 
model in the dichotomous Markovian field is possible by analogy with the 
consideration pursued in Sec. 15.41 The result is formally the same as in 

19 Reminder: we consider the case A r = 0, or fcgT ^> HA r . 
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Eq. with f (t - f) and g (t - t') given but in Eq. (jTUl) (with A(t) 

= const) [101]. An interesting feature is that for e = the equations for 
the average (cr z (t)) and the correlator (a(t)a z (t)) are decoupled. Moreover, 
in the dissipation-free case, Q(t) = 0, the solution of equation for (cr z (t)) 
with the initial condition (<r z (0)) = 1 yields the same result as in Eq. fl40(l 
with the following substitutions implemented, i.e., (a z ) — > (<J X ), o — > A, 
A — > eo- This finding provides a rather nontrivial cross-check of the validity 
of different methods of stochastic averaging. 



6.3.3 Electron transfer in fast oscillating periodic fields 

Let us next focus on the case with strong and fast periodic driving fields 
e(t) = Acos(Qt + <fo), yielding 

oo . 
n=—oc 

(W 21 (e )) e = -A 2 J n(^y kBT FC(e -nQ). (112) 

n=— oo 

This result of the Golden Rule type for the nonadiabatic electron transfer 
(ET) rates in strong periodic fields has been derived in [83] and independently 
in [89]. In particular, the quasi-static (Gaussian) approximation for FC{uj) 
for k B T > fku c with K(t) replaced by K(0) w 2k B TX/h 2 in Eq. (fT05jl leads 
in absence of driving, i.e. A=0, independent of the detailed structure of J{u), 
to the celebrated Marcus-Dogonadze-Levich rate expression [184, 185] for the 
ET rates with 

This approximation is suitable for a thermal bath with a low frequency cut- 
off and in the high-temperature limit, e.g., for polar solvents. This presents 
a semiclassical limit for the Franck-Condon factor. If some high-frequency 
(quantum) vibrational mode ujq couples to ET with the coupling constant 
kq in addition to the low-frequency vibrations, being of relevance for ET in 
molecular aggregates, then a different model for FC(uo) is more appropriate, 
namely [8, 186], 



Pj — phiJQ 



p=—oo 

X 6XP ( 4A^T ) (U4) 
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where S = (^) 2 , D = Scoth(^), x = 5/sinh(^), and I p (x) is the 
modified Bessel function. The periodic driving may induce an inversion of 
ET transfer direction and modulate ET transfer rates by orders of magni- 
tude. This has been theoretically predicted in [83,89] for both of the above- 
mentioned models of FC(u). 20 

6.3.4 Dichotomously fluctuating tunnelling barrier 

Another relevant situation involves the case of fluctuating tunnelling ma- 
trix element A(t) and constant energy bias e(t) = eo = const. In the su- 
perexchange picture of ET this corresponds to a physical situation where 
the stochastic dynamics of the bridge states, which mediate the electron 
transfer between the donor and acceptor molecules, introduces an explicit, 
stochastic time-dependence into A(t). Generically, this corresponds to a fluc- 
tuating tunnelling barrier. In the case of dichotomous Markovian fluctuations 
A(t) = Ao+Aa(t), the stochastic averaging of the NIBA master equation can 
be done exactly [67]. Towards this goal one makes use of the Shapiro-Loginov 
theorem (|%3~j) and the following exact decoupling property [152,153,187]: 

(a(t)a(t>)o z (t>)) = (a(t)a(t'))(a z (t')) . (115) 

Applying these two theorems and using the DMP property, a 2 (t) = 1, the 
averaging of the GME yields the following exact results [67]: 

J' ([A 2 + AV^]/(t-t')(^(t')) 
A A[1 + e-^]f(t - t')(a(t')a z (t')) 
[A 2 + AV^^Mt - t')) dt', (116) 

u(a(t)a z (t)) - jT ([A 2 + A 2 e^-*')]/(t - t')(a(t')a z (t')) 

A A[1 + e-^]{f(t - t')(a z (t')) + g(t - t')})dt' , 

= exp[— Re Q(t)] cos[lmQ(t)} cosfeot], 

= exp[-ReQ(t)]sin[ImQ(t)]sin[e t]. (117) 

20 The use of an improved perturbation theory in A in the case of fast fluctuating fields 
does not imply that the Golden Rule rates cannot be enhanced by such nonequilibrium 
fields. A large enhancement of the forward (backward) rate can occur, e.g., when the ab- 
sorption of n photons helps to overcome the corresponding forward (backward) activation 
barrier of the thermally-assisted incoherent tunnelling. For example, for the generalised 
Marcus rates a condition is cq X/h± nhVl = 0, with the field amplitude A chosen such 
that the probability J% (A/£l) of the corresponding reaction channel is maximised. 



where 



+ 
+ 



/(*) 
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For the case of vanishing dissipation, Q(t) = 0, and for A = 0, the solution of 
this integro-differential equation for (a z (t)) for the initial condition (er 2 (0)) = 
1 yields the same result as in Eq. (}35|) . This agreement provides an additional 
test for the mutual consistency of different methods of stochastic averaging 
used here. 

Furthermore, in the absence of dissipation the rate of incoherent relax- 
ation exhibits a resonance-like feature as a function of the frequency v of the 
barrier fluctuations. Namely, a resonance occurs when v matches the tran- 
sition frequency eo, i.e. v = e (see Eq. fp?7|) in Sec. 12.3.2)1 . This presents a 
physical stochastic resonance, which should not to be identified with a well- 
known phenomenon of noise-assisted Stochastic Resonance [188]. It occurs 
when a stochastic frequency of the driving matches an eigenfrequency of a 
quantum transition. In the presence of dissipation, this resonance feature 
is maintained, but becomes modified. Namely, the resonance can occur at 
v = \eo±X/h\, rather than at v = eo [67]. This resonance is responsible for the 
interesting phenomenon of a stochastic acceleration of dissipative quantum 
tunnelling which is predicted by the theory [67]: For the case that A := A , 
when A(t) fluctuates between zero and 2A the rate of incoherent trans- 
fer can exceed that for the static tunnelling barrier with tunnelling coupling 
strength A(t) = 2Ao =const. At the first sight, this effect seems paradoxical; 
it must be remembered, however, that the considered noise is nonequilibrium 
and it is capable of pumping energy into the system enhancing thereby the 
rate of incoherent quantum tunnelling. Unfortunately, for the parameters 
typical for molecular ET the experimental conditions for this effect to occur 
can barely be met experimentally because the required frequency v is too 
high. Nevertheless, this fact does not invalidate the principal possibility of 
the discussed effect for some other physical systems in view of the generality 
of the model set-up. 

In contrast, when A(t) fluctuates very slowly on the time-scale of decay of 
kernels f(t) and g(t), which corresponds roughly to the inverse of the width of 
corresponding Franck-Condon factor FC(u), then our theory predicts - after 
use of the Markovian approximation - the known results which corresponds 
to the approximation of a dichotomously fluctuating rate [189], see also dis- 
cussion in Sec. 15.31 The corresponding problem of such fluctuating rates is 
known under the label of dynamical disorder and can be met in quite differ- 
ent areas of physics and chemistry [190] . Depending on the relation between 
the stochastic frequency v and the values of transfer rates corresponding to 
the "frozen" instant realizations of A(t), the transfer kinetics can exhibit 
different regimes of a (i) quasi-static disorder, (ii) an averaged rate descrip- 
tion, and (iii) a gated regime [67]. In the latter case, the mean transfer time 
becomes locked to the autocorrelation time of the fluctuations [68,191]. 
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The influence of strong laser fields on the electron transfer with nonequi- 
librium dynamical disorder [86], or driven by nonequilibrium conformational 
fluctuations [68] has been studied within the obtained NIBA master equa- 
tion approach in Refs. [68,86]. In particular, it has been shown there, that a 
strong periodic field can induce a turnover between the nonadiabatic regime 
of electron transfer and a gated regime. Moreover, the direction of elec- 
tron transfer in the gated regime can be inverted, whereas the mean transfer 
time remains chiefly controlled by the nonequilibrium stochastic fluctuations 
and it is not influenced by periodic field [86]. These theoretical predictions 
discussed here are still awaiting their experimental realization. The area of 
chemically gated, or chemically driven electron transfer [192], that is the elec- 
tron transfer controlled by nonequilibrium fluctuations due to spontaneous 
release of energy by breaking some energy-rich chemical bonds (e.g., due to 
the ATP hydrolysis), is currently still in its infancy [192]. 

7 Quantum transport in dissipative tight- 
binding models subjected to strong 
external fields 

A salient application of our theoretical efforts relates to charge and parti- 
cle transfer in spatially extended molecular structures. It can be described 
within a model similar to the Holstein model of a molecular crystal [167, 193]. 
Namely, one considers a molecular chain using the assumption that only 
one energy, namely a lowest, unoccupied molecular orbital(LUMO)-level de- 
scribes the electron transfer (or highest occupied, HOMO-level in case of a 
hole transport) per molecule, or molecular group. These energy levels are 
coupled to the local intramolecular vibrations which are thermalised. The 
transferring particle is delocalised due to a tunnelling coupling between the 
nearest neighbours. The intersite coupling between the intramolecular vibra- 
tions is however neglected (like in the Einstein model of optical phonons), i.e. 
the electron (or hole) energy levels in neighbouring molecules (or molecular 
groups) are assumed to fluctuate independently. In other words, one assumes 
uncorrelated identical thermal baths formed by vibrational degrees of free- 
dom of each molecule in a molecular chain. Such a model is close in spirit to 
one used for exciton transfer within a stochastic Langevin description [40]. 
In the approximations employed below, this model becomes equivalent to the 
model of a Quantum Brownian Motion within a single band, tight-binding 
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description. In an external electric field S(t), the latter one reads [3,4]: 



hA x 



H TB (t) = -— J2 (\n)(n + l\ + \n+l)(n\)-e£(t)x + H m , (118) 



2 

n=—oo 



2 z — ' Vvtii V m,;a>, 



| g; - — 2^ 



where x = a n\n) (n\ is the operator of the coordinate (within the single 
band description). The model in Eq. ()118j) can be derived from a different 
perspective than the Holstein model, namely, by starting out from a model of 
Quantum Brownian Motion in a periodic potential [2,3,80] and by restricting 
the corresponding consideration to the lowest band for the tunnelling particle 
in the deep quantum regime. We consider this model in the limit of a strong 
coupling by applying the small polaron transformation which now reads U = 
exp[— ix P/h], P = J2 i Cipi I (rriiUjf) . In the polaron basis, the Hamiltonian 
reads 

hA °° 

H TB (t) = (\n)(n + l\ + \n + l)(n\)-e£(t)x 



2 

n=— oo 



2 



- E (i\n}(n + l\ + h.c.) + -J2[^: + ^Q^ (H9) 

n=— oo i 

where A r = A(e~ mP / h )B = Ae^ 0,2 ^ 2 ^ 13 ^ 2 ^ is the renormalised tunnelling 
coupling (polaron band width), Qi := UqiU~ l = ^x are displaced bath 

oscillators and ^ = MAe~ iaP ' h — A r \ is the quantum random force operator in 
the polaron basis which is considered further as a small perturbation. Note 
that x is not changed. Assuming a strong Ohmic dissipation with a > 1 yields 
A r = at T = OK and for £ (t) = 0. This indicates the celebrated localisation 
phase transition [176,177], which alternatively can also be interpreted as a 
polaron band collapse. In the presence of a constant electric field and/or for 
T > this localisation transition occurs for smaller values of a. Given our 
case of strong coupling, the transport occurs predominantly via incoherent 
tunnelling hops between the nearest sites of localisation. As a side remark, 
we note that also in the dissipation-free case the Bloch band can collapse 
in presence of strong periodic fields [194], known as the effect of dynamical 
localisation [195]. Use of Eq. (jHTjl for the case in (|119j) with A r = yields 
for the diagonal elements of the reduced density matrix a set of coupled 
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generalised master equations 



Pnn(t) = / {W (+ \t,r)p n _ ln _ 1 (r) + W { -\t,r)p n+ln+1 (T) 
Jo 

- [W^(t,r) + W^\t,T)]p nn (r)}dT (120) 
with memory kernels 

W {± \t,r) = ^A 2 e' RcQ( '- r) cos[Im Q(t - r) T ^ / (121) 

The very same equations are obtained in the NIBA approximation of the 
path-integral approach [196]. The Holstein like model which has been dis- 
cussed at the beginning of this Section yields in similar approximations the 
same set of GME's (with a trivial renormalisation of the coupling constant 
in the identical bath spectral densities J n {ui) = J{u))) [69]. The stationary 
electrical current carried by one particle reads j = elim^oo -^(x(t)), where 
(x(t)) = ciJ2n n Pnnit) denotes the mean particle position in the considered 
infinite chain. It obeys (this result follows immediately from Eq. ()120j0 



, / [W + (t,r) - W-(t,T)]dr. (122) 
Jo 



This current in ()122|) still needs to be averaged of the stochastic field re- 
alizations. 21 This objective is again reduced to the averaging of an effec- 
tive Kubo oscillator which can be done exactly for many different models of 
stochastic driving. We decompose the electric field £{t) into the sum of the 
mean, or constant field So and a fluctuating, unbiased component £{t), i.e. 
£{t) = £q + £{t). The resulting expression for the averaged current j(£o) can 
be put into two equivalent forms. First, it can be written in terms of a time 
integral [70,102], 

POO 

j(£ ) = eaA 2 / exp[-Re Q{r)\ sin[Im Q(T)}lm[e iea£ ° T/h (S(r))]dT, (123) 
Jo 

where (S(t)) is given in Eq. ()108|) with e(t) = ea£(t)/h and Q(t) in Eq. 
()105j) . Alternatively, the current expression can be given as a frequency 
convolution in a spectral representation form, i.e., 

POO 

j{£ ) = / j dc {u)I(ea£ /h- u)du, (124) 



2 In the case of periodic driving, this additional averaging is obsolete by defining the 
current in a self-averaged manner as j — e linif^oo \ x Wi 
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where 



j dc {u) = -eaA 2 (l - e-W")FC{u) . (125) 

I{ui) denotes the spectral line shape corresponding to (S(r)). The dc- 
current obeys the symmetry property jdc(— u>) = —jdc(u) which is imposed 
by the thermal detailed balance symmetry, FC(— uj) = e~ hl3uJ FC(uj) with 
FC{uj) in (jllOj) . It is important to note that the averaged current in (jl24j) 
does not obey such a symmetry requirement. 

7.1 Noise- induced absolute negative mobility 

As a first application of the above results we consider the phenomenon of 
ANM, or absolute negative mobility, where the transferring particles move 
around zero bias in opposite direction to the average applied force. This 
effect has been anticipated for semiconductors in strong periodic fields almost 
thirty years ago using a Boltzmann equation approach [197,198]. The first 
experimental realization was obtained in 1995 for semiconductor superlattices 
[199]. The corresponding experimental results were seemingly consistent [199] 
with a mechanism of incoherent sequential tunnelling like one just described. 
The occurrence of the ANM phenomenon for a sinusoidal driving within the 
considered dissipative tight-binding model has been demonstrated in Ref. 
[196]. 

The question we addressed in Ref. [69] within a Holstein like model was 
whether an external stochastic field can also induce ANM. The occurrence of 
such noise-induced ANM has been shown for dichotomous Markovian fields. 
ANM presents a multi-state analogy of the effect of inversion of populations 
in TLS described in Sec. El A simple criterion for ANM to occur can be 
found within the quasi-static approximation for the spectral line shape 
For a symmetric dichotomous field £(t) = (ho~/ea)a(t) with the inverse au- 
tocorrelation time v, this quasi-static approximation holds whenever a ^> u, 
being almost always the case in the relevant regime of parameters even if 
the field fluctuations are fast on the time-scale of the charge transfer. Then, 
« \[8{u-a)+8{u-a)] and j(£ ) = ^[j dc (ea£ /h-a)+j dc (ea£ /h+a)}. 

Given the symmetry property, jdc(~ cr ) — — Jdc( cr ), one can conclude that 
the phenomenon of ANM will occur in any such system with the static 
cur rent- volt age characteristics idc{.o~) assuming a maximum at some a max 
which is complemented by a corresponding regime of differential negative 
conductance occurring for a > a max . Then, j(£o) < for a sufficiently 
small static force, e£ > 0, whenever a > o max [69], i.e. whenever the 
charge transfer is driven into the regime of negative differential conductance 
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by some appropriately chosen alternating, two-state stochastic fields. This 
mechanism is quite general and robust. It does not depend on the details 
of the dissipation mechanism. In particular, for the Gaussian FC(uj) in Eq. 
(jllHj) . we obtain 



j dc (ea£ /h) 



it eaA 2 h 



2 y/irXksT 



cxp 



A 2 + (eaS 
AXk B T 



sinh 



eaSo 
2k B T\ 



(126) 



This corresponds to a (nonadiabatic) small polaron conductance [167, 200] 
with the differential mobility, /i(£o) — , obeying in the linear response 
regime 



Mo) 



7T 



ea 2 V 2 



2W P h{k B Ty/ 2 



-W p /2k B T 



(127) 



where W p = A/2 is the polaron binding energy and V = HA/2. For this 
nonadiabatic small polaron model the regime of negative differential mo- 
bility occurs for £ > S max with E max defined implicitly by the equation 
eaEmax = 2W P coth( g ffgfr ) . Quasi-one-dimensional systems exhibiting this 
small polaron conductance (in the nonadiabatic ET regime with respect to 
A) can be considered along with the semiconductor superlattices as possible 
candidates to exhibit the phenomenon of noise-induced ANM experimentally. 
A finite photo-induced small polaron mobility of the hole type is found, for 
example, in columnar liquid crystals [201-203]. We then estimate the value 
of S max for these systems with the lattice period of about a = 0.35 nm to be 
in the range of 5 • 10 6 V/cm, which is rather large. For superlattices with a 
larger period a, S max can be much less [199]. Basically, this crucial quantity 
is determined by two factors: (i) the width of FC{uj) due to multi-phonon 
transitions (it depends on the precise mechanism of dissipation and should 
be made as small as possible) and (ii) the lattice period a (it should be en- 
gineered as large as possible). These criteria can serve as a useful guides in 
identifying the appropriate experimental materials. 



7.2 Dissipative quantum rectifiers 

Yet another intriguing application is provided by the fluctuation-induced 
quantum transport in the absence of a mean electric field, Sq = 0. Similar 
nonequilibrium phenomena are known under the notion of Brownian motors, 
or Brownian ratchets [71-81]. The first case of a quantum ratchet in a pe- 
riodic spatially asymmetric (ratchet)-potential was studied theoretically in 
Ref. [204] within a semi-classical approach and for an adiabatically varying 
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driving field. In Ref. [70,102,104], we put forward periodic dissipative nona- 
diabatic quantum rectifiers [209] operating in the absence of spatial asym- 
metry. The current is produced by a nonlinear transport mechanism due 
to an interplay between equilibrium quantum fluctuations and an unbiased, 
but asymmetric nonequilibrium external noise [70]. Likewise, an asymmet- 
ric periodic driving of the harmonic mixing type can be used instead of the 
nonequilibrium noise [102, 104]. Our rectifier behaves genuinely quantum 
mechanically and corresponds to the case of a strong dissipation when the 
transport mechanism is incoherent and the transport proceeds by incoherent 
tunnelling hopping as outlined above. The origin of the resulting current 
can be traced to Eq. (ITIZHJ) and Eq. Namely, j(0) ^ 0, when (S(t)) 

assumes a complex values, i.e., Im(5'(r)) ^ 0. This corresponds to a comple- 
mentary criterion which follows from Eq. ([124)1 . Namely, j(0) 7^ 0, when the 
corresponding spectral line is asymmetric, I{—oj) 7^ 

In particular, this is the case of asymmetric dichotomous field of zero 
mean, cf. Sec. 12.3.11 and Eq. (J2HJ), which takes on the (frequency scaled) 
two discrete values ea£ 12 /h = ei j2 = =Fcre Tfe//2 , where b characterises the field 
asymmetry and a is the (scaled) rms of field fluctuations. The emergence of a 
finite current in this case can readily be seen in the quasi-static approximation 
of I(u>) for a 3> u, I{u) « Pi5{u — aer h / 2 ) + p 2 5(uj + oe h l 2 ) with p lj2 = 
|£2,i|/(|ei| + £2)- In this adiabatic (with respect to driving) approximation, 

j(0) = P2Hc{ve b ' 2 ) - Pl ] dc {ae~ b ' 2 ). (128) 

In the semiclassical high-temperature approximation for FC (u) in Eq. ()113j) 
(this corresponds to noise-driven small polaron transport), one can see that 
the rectification current appears as a nonlinear response to the external, 
unbiased on average driving. Namely, to the lowest order, the current is 
proportional to (£ 3 (t)), j(0) oc (£ 3 (t)) « ba 3 {b <^ 1), with a nontrivial 
prefactor. Moreover, the current flows into the direction of (e 3 £ 3 (t)), which 
is the direction of the larger force realization, if the applied random force 
is sufficiently small. With an increase of the noise rms a the current can 
however reverse its direction. In the considered approximations and for a 
small driving asymmetry 6 < 1 this occurs when a exceeds some maximum 
associated with FC(u). Thus, this change of the current direction from 
the expected to physically counter-intuitive direction is closely related to 
the mechanism of noise-induced absolute negative mobility, as it has been 
outlined in the previous subsection. 

Moreover, the current can flow in the physically counter-intuitively direc- 
tion also for small applied forces when the coupling strength a is sufficiently 
small. For T = a very insightful approximate analytical expression can 
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analytical result — 
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Figure 7: Noise induced rectification of current for an asymmetric dichoto- 
mous driving field vs. rms of field fluctuations at different temperatures. The 
set of used parameters is indicated in the figure. 

be obtained in the adiabatic limit for driving and in the lowest order of the 
asymmetry parameter b. Namely, assuming an Ohmic friction mechanism 
for the thermal bath with an exponential cutoff u c , FC(uo) can be exactly 
evaluated at T = to yield FC(uo) = u J [2a) (^) 2a ~ 1 exp[-uj/uj c ]Q(uj), where 
Q(u>) is the Heaviside step function [2]. For b <C 1 in ()128j) . this then yields 
the averaged, zero bias current value, reading for a > 1/2: 

M ~ 6- ^ A -T'U -l-—) exp(-aM) . (129) 
JK ' 2 uj c T (2a) \ujJ V 2uJ PV 1 ! K J 

The result in Eq. ()129|) predicts that for a < 1 the current flows into 
the physically counter-intuitive direction. Furthermore, for a > 1, the recti- 
fication current flows first in the expected, natural direction, but it changes 
subsequently its direction for a > a* = 2 (a — l)u c . Moreover, the absolute 
value of current has two maxima at <J max = (2a — 1 ± \j2a — 1)uj c for a > 1 
and one maximum at a max = (2a — 1 + \/2a — V)uj c for 1/2 < a < 1. Fur- 
thermore, the current diminishes for large a. In the low-temperature limit, 
all these features are in the remarkable agreement with the numerical eval- 
uation of Eq. ()12Hj) in Ref. [70]. A related comparison is provided in Fig. 
[7|for a = 2. For k^T = 0.0lhu) c , the agreement is indeed excellent, except 
for very small values ha &bT. On the scale of a variation used in Fig. 
[7| the rectification tunnel current seems be maximal at T = for most a. 
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From this point of view, the rectification results from an interplay between 
the zero-point quantum fluctuations and the nonequilibrium noise, i.e. has a 
manifestly quantum origin. 

For larger values of a than those depicted in Fig. Q the rectification 
current can be enhanced for sufficiently small <j/uj c by increasing temperature 
and go through a maximum, exhibiting thereby the phenomenon of Quantum 
Stochastic Resonance [205-208] in the nonlinear current response [70]. 

7.3 Limit of vanishing dissipation 

In the limit of vanishing dissipation Q(t) — ► 0, the result in Eq. ()123j) 
predicts that j(0) = 0, independently of the form and strength of driving. 
This prediction should be considered, however, with care since the result in 
Eq. ()123|) is not valid for very small a and T, because we have assumed 
throughout incoherent transport regime where either A r = 0, or the tem- 
perature is sufficiently high, ksT ^> HA r . Nevertheless, the dissipationless 
single-band, infinite tight-binding model can be solved exactly in arbitrary 
time-dependent fields [70,103,104,195,210,211]. The corresponding exact 
solution for the current then shows [70, 103, 104] that the stationary current 
is forced to vanish identically by the stochastic fluctuations of driving. Put 
differently, in the absence of quantum dissipation such a rectified current can 
exist at most as a transient phenomenon. As a matter of fact, the stationary 
rectification current within the single-band tight binding description is due 
to a nonlinear interplay of quantum dissipation and external nonequilibrium 
forces. Its origin presents a highly nonlinear and nonequilibrium statisti- 
cal effect. This result does not hold for more general situations. For the 
case of the full potential problem, with its intrinsic interband transitions, 
a finite, stationary current can be generated even in the absence of dissi- 
pation; it results as a dynamical effect due to an interplay of a nonlinear 
dynamics and the breaking of some space-time symmetries by the driving 
mechanism [103,212,213]. 

It must be emphasised, however, that the full potential problem has little 
relation to the electron transport in molecular chains which is our main 
focus here. This is so, because the tight-binding description emerges for the 
electron (or hole) transport processes in molecular systems in a very different 
way, being not the result of a truncation of a full potential problem to the 
description within the lowest band only. 
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7.4 Case of harmonic mixing drive 

Another instance of quantum rectifiers in presence of dissipation is realized 
with a harmonic mixing driving [102,214,215], 

£{t) = E 1 cos(fit) + E 2 cos(2ftt + 0), (130) 

with the driving strengths Ei,E 2 , angular frequency Q and a relative phase 
0, respectively. This model seems more promising and readily can be imple- 
mented with experimental realizations. 

The corresponding expression for (S(t)) reads [102] 

oo 

(S(r))= ^ fc (26sin(fir/2))j fc (e 2 sin(fir))e-^ +7r / 2 ), (131) 

k=— oo 

where £i t 2 = eaEi t2 /(hQ) and J n (z) are standard Bessel functions. With its 
help the current in Eq. ()123j) can be evaluated numerically for the Ohmic 
model with the exponential cutoff, where the exact analytical expression for 
Q(t) is available [2,3,69]. Independently of other parameters the current van- 
ishes identically for = 7r/2,37r/2, where lm(S{r)) = exactly. Otherwise, 
the current can be different from zero. For sufficiently high temperatures and 
weak fields applied, j(0) oc (£ 3 (t)) = jEf E 2 cos(0) with a nontrivial quan- 
tum prefactor. At T = 0, the current response is not analytical in the driving 
amplitude. Unfortunately, in this case we do not find a simple approximate 
analytical expression for the current like the one in Eq. (|129jl . Some nu- 
merical calculations [102], see also in Fig. El reveal a series of nontrivial 
features as the occurrence of current inversion and the occurrence of current 
maxima similar to the the case of stochastic dichotomous driving. More- 
over, in the case of harmonic mixing driving the direction of the rectification 
current can be conveniently controlled by the phase 0. For a sufficiently 
large dissipation strength a, the rectification current response can also ex- 
hibit a Quantum Stochastic Resonance feature [188,205-208,216-218], i.e., 
it exhibits a maximum versus the temperature T. An experimental realiza- 
tion of the dissipative quantum rectifiers in the studied incoherent tunnelling 
regime can be expected for the semiconductor superlattices [219] and for a 
small polaron like transport in molecular chains. 

8 Summary 

We have surveyed, extended and justified in great detail the results of recent 
research which relates to the quantum dynamics with fluctuating parameters. 
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Figure 8: The rectified quantum current induced by a harmonic mixing drive 
is depicted vs. the strength of the first harmonic at different values of the 
coupling strength a. The strength of the second harmonic is held fixed. The 
used parameters are given in the figure. 



The nature of those fluctuations, that usually stem from the influence of 
externally applied fields, or intrinsic degrees of freedom, is either of stochastic 
or coherent deterministic origin. Our general findings are applied to specific 
situations, encompassing solid state applications, but mainly we did focus on 
the theme of driven electron transfer in molecular aggregates. 

A first major result is the derivation of the generalised non-Markovian 
master equations and the generalised Redfield equations for quantum sys- 
tems composed of a finite, discrete number of states which are subject to 
the influence of external either stochastic, or periodic fields. The result- 
ing kinetic equations allow one to study a rich variety of different physical 
problems within a unified framework. In simplest cases, the relevant part 
of the reduced dynamics is either described by the Markovian balance equa- 
tion of the Pauli master equation type [220] which is generalised to include 
the influence of explicit time-dependent, external field manipulations into the 
quantum transition rates; those becoming therefore functionals of the driving 
field, or by its generalisation which accounts as well for the memory effects 
in the corresponding dissipative kernels. 

In the case of fast (on the time scale of the averaged relaxation process) 
fluctuating, or oscillating fields these quantum kinetic equations can be aver- 
aged. The relaxation transitions can be described by the averaged quantum 
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transition rates of the Golden-Rule type. These averaged transition rates, 
however, generally do not satisfy the detailed balance condition at the tem- 
perature of the thermal bath. This violation of the thermal detailed balance, 
being induced by the nonequilibrium driving fields, in turn paves a roadway 
for identifying several intriguing nonequilibrium phenomena. 

An important case is provided by symmetric dichotomous driving fields 
for which our approximate theories and considerations can be made rigorous, 
tested and reaffirmed [67, 100, 101] because the corresponding averaging can 
be performed exactly. 

The problem of averaging a quantum dynamics in stochastic fields mod- 
elled by non-Markovian processes of the continuous time random walk type 
with a discrete number of states (with Markovian processes emerging as a 
particular limiting case) has been investigated in the absence of dissipation in 
Sec. 2. Using a classical stochastic path integral approach, we obtained some 
general exact results on the averaging of quantum propagator of the driven 
quantum system over the stationary realizations of such non-Markovian jump 
processes. In particular, the exact result for the Laplace transform of the cor- 
respondingly averaged quantum propagator has been obtained. This novel 
result bears the potential for future applications since it opens a way for a 
rigorous study of an extreme case of l/f a noise, implying long-range tem- 
poral correlations, where standard perturbation theory is expected to fail. 
As a first important application, we obtained the spectral line shape of the 
corresponding Kubo oscillator and the Laplace-transformed averaged evolu- 
tion of a spin 1/2- system that is driven by a symmetric alternating renewal 
process possessing an arbitrary distribution of the residence times. This im- 
plies a very broad class of autocorrelation functions including those which 
correspond to noise sources with 1/ f a power spectrum. This general result 
is shown to reproduce the known solution in the Markovian limit, i.e. when 
the the residence time intervals are exponentially distributed. 

Starting out with section 4, we have investigated the combined influence of 
fluctuating parameters and dissipation on the evolution of the corresponding 
quantum dynamics. In doing so, we considered both the role of fluctuations of 
the energy bias, and/or the effects of a fluctuating intersite tunnelling matrix 
element. A generalised master equation was obtained which corresponds to 
the known NIBA approximation obtained within the quantum path-integral 
approach. The obtained master equation was averaged exactly both over 
dichotomous fluctuations of the energy bias and dichotomous fluctuations of 
the tunnelling coupling, i.e. the case with a fluctuating tunnelling barrier. 
These results have been used to study a rich repertoire of nonequilibrium 
phenomena for electron transfer in condensed media with dynamical disorder 
and possibly being driven by stochastic or deterministic, coherent periodic 



67 



laser fields. 

In section 7 we have studied the quantum transport in extended quan- 
tum systems within a tight-binding description, with the dynamics being 
subjected to a strong system-bath coupling and weak tunnelling, i.e. in the 
limit of an incoherent hopping regime. A general result for the quantum- 
noise assisted transport current, being averaged over the field fluctuations, 
has been derived. The corresponding expression is shown to be equivalent 
to the NIBA approximation result of a corresponding quantum path integral 
treatment. 

Our theory for dissipative systems with fluctuating parameters predicts 
scores of interesting nonequilibrium phenomena that are the result of a stun- 
ning interplay between equilibrium quantum fluctuations and nonequilibrium 
perturbations. A few noteworthy such effects are: (i) the suppression, or ac- 
celeration of quantum transition rates by many orders of magnitude; (ii) 
a noise-induced enhancement of the thermally assisted quantum tunnelling; 
(iii) the inversion of populations in the spin-boson model; or (iv) a noise- 
induced absolute negative mobility in quantum transport. We further elab- 
orated on the theme of dissipative quantum rectifiers. Several of these novel 
predictions are presently being investigated in a number of research groups, 
both theoretically and experimentally. Our research in particular also im- 
pacts such timely activities like the investigation of the electronic transport 
in infrared laser driven molecular wires [221]. Here, the fermionic thermal 
baths are provided by the electronic reservoirs in the leads and the elec- 
tron transport through the wire is mainly coherent. This corresponds to 
the regime of a weak dissipation within our approach, being opposite to the 
regime of incoherent tunnelling. The role of the size, inter-electrode coupling 
effects, etc. [222] as well as the inelastic Coulomb repulsion effects [223] , are 
also important for molecular wires. This brings about further complications 
that still need to be investigated theoretically with greater detail. Experi- 
mental progress is presently also forthcoming [224]: This particularly holds 
true for quantum Brownian motors and quantum rectifiers as witnessed by 
the exemplary set of recent experimental studies [225-230]. We share the 
confident belief that this research topic will remain flourishing and, more- 
over, will invigorate the readers in pursuing their own future research in this 
area. 
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